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ABSTRACT 

We derive the extinction curve towards the Galactic Center from 1 to 19 /im. We use hydrogen 
emission hues of the minispiral observed by ISO-SWS and SINFONI. The extinction free flux reference 
is the 2 cm continuum emission observed by the VLA. Towards the inner 14" x 20" we find an extinction 
of A2. 166pm = 2.62 ±0.11, with a power- law slope of a = —2.11 ±0.06 shortward of 2.8 /zm, consistent 
with the average near infrared slope from the recent literature. At longer wavelengths, however, 
we find that the extinction is grayer than shortward of 2.8 /im. We find it is not possible to fit 
the observed extinction curve with a dust model consisting of pure carbonaceous and silicate grains 
only, and the addition of composite particles, including ices, is needed to explain the observations. 
Combining a distance dependent extinction with our distance independent extinction we derive the 
distance to the GC to be Ro = 7.94±0.65 kpc. Towards Sgr A* (r < 0.5") we obtain Ah = 4.21 ±0.10, 
Aks = 2.42 ± 0.10 and Al' = 1.09 ± 0.13. 

Subject headings: Galaxy: center - Galaxy: fundamental parameters- ISM extinction, dust 



1. INTRODUCTION 

Knowledge of the e xtinction in the infrared 
(ISchultz fc Wiemeri [19751 : ICardelli et aP [19891 : IMathisI 
|1990( ) is important for obtaining the intrinsic lumi- 
nosities of highly extincted objects. In addition, the 
extinction curve provides important c onstraints on the 
prope rties of interstellar dust (see e.g.. iCompiegne et al.l 
1201 1|) . Using the extinction curve together with the 
dust emission it is possible to derive the composition 
and sizes of interstellar dust grains. In this way it has 
been found that interstellar dust is mai nly compose d 
of silicate and carbonaceous dust grains (IDraind 120031) ■ 
However, it is still uncertain if composite particles 
are also required jM athis 1996; Li fc Gr cenbcrg 1997; 
iWeingartner fc DraineH200ll : IMAo et al.ll2004l ). Includ- 
ing, for example ices and voids in addition to the basic 
dust grains. 

The Galactic Center (GC), invisible in the optical, 
is one well-known example o f a highly extincted region 
([Becklin h NeugebaueH[l968l ). Because of this, it should 
be relatively easy to measure the infrared extinction, 
without the large relative errors impl icated in using re- 
gions of small absolute extinction. iRieke fc Lebofskvl 
}r985), for example, used infrared photometry of a few 
red supergiants present in the GC together with a vis- 
ible to near infrared extinction law obtained outside of 
the Galactic Center, in order to derive a universal in- 
frared extinction law (Figured]). This measurement was 
later slightly imp roved by addi ng NICMOS near infrared 
measurements bv [Rleka jTUP). 

Another method of measuring the IR extinction is 
to u se a stellar popu l ation o f known infrared luminos- 
ity ([Nishivama et al.l l2006bl ). such as the red clump 
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for which the absolute luminosity is measured locally 
(|Groenewegen| [2008l . For this method, however, it is also 
necessary to know a precise distance to the object, which 
makes the method difficult for most objects. In the case 
of the GC the dista nce is however well known (Rcid 19931 
iGenzel et al.l 120101 ) . Using this method Schodcl e t ahl 
(12010D obtained A^ = 4.35 ± 0.18, A^^ = 2.46 ± 0.12 
and Kl' = 1.23 ± 0.20 towards Sgr A*. 

As an alternative to methods involving stars, one can 
also use nebular hydrogen lines to measure extinction. As 
a case in point. .Lutz et al. (199(1) and[Lutz (1999) used 
ISO-SW S spectra of the minisp iral (a bright HII region in 
the GC; iLo fc ClaussenI [T983i i to derive extinction from 
2.6 to 19 /xm. Typically this method is more precise, since 
the intrins ic uncertainty of the rela tive line strengths in 
this case ([Hummer fc Storevlll987D is smaller than the 
spectral uncertainty of stellar emission. An additional 
advantage is that the lines are much narrower than the 
bandpasses of the broadband filters, so that the effective 
wavelength is known a priori, while for broadband filters 
it is extinction dependent. Furthermore, there are many 
more lines in the IR than the number of broadband filters 
in frequent use, such that it is easier to obtain a well- 
sampled extinction curve, that includes also extinction 
features. 

Although it is only possible, using infrared lines alone, 
to obtain relative extinction measurements within the 
infrared, stellar methods often have the same problem. 
Absolute extinction calibration is made possible by com- 
paring with longer wavelengths where the extinction is 
negligible. The absolute extinction towards the GC for 
th e Paschen-g line at 1.87 ^m was obtained in this way 
m IScoville et all ([200l using 6 cm continuum data. 

The extinction towards Sgr A* from the various mea- 
surements mentioned above is shown in Figure [1] How- 
ever, there is apparently some discrepancy between the 
different studies. The uncertainty in extinction creates, 
in turn, uncertainty in determining the intrinsic lumi- 
nosit i es of young mass ive stars in the GC ([Martins et al.l 
[2007t[fVitz et al.|[2010l) . as well as the intrinsic luminos- 
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Fig. 1. — Extinction towards Sgr A* from the literature; triangles 
represent values derived from stars, boxes values derived from gas 
lines. Most values were obtained using slightly different regions 
around Sgr A*. For the comparison, we convert the values to 
the direct sight line towards Sgr A* (i. e. the centra l 0.5" , see 
Section |5.3|I using the extinction map oflSchodel et al.I II 2010I ). In 
the case of Rickc & Lebofskv (1985) we use the extinctio n law of 
[Rieke fc Lebofskvi (.1985. ) and the absolute extinction of iRieke et al.l 
H1989I 1. 



ity of Sgr A* (iGenzel et al.l 120031 ) . In addition, there is 
no study that covers the full wavelength range, and only 
a few studies of this kind are available outside of the 
Galactic Center. 

In this paper we use the emission of the niinispiral from 
1.28 to 19 /im to derive relative extinction values from 
line emission, which are compared to 2 cm radio contin- 
uum data to compute the absolute extinction. In Section 
2 we present our data, and discuss the extraction of line 
maps and flux calibration in Section 3. In Section 4 we 
derive the extinction. We discuss the results in Section 
5. We compare the results with other works for the GC 
and for other parts of the Galaxy and with theoretical 
dust models. We also use the data for estimating the 
distance to the Galactic Center, Rq. We summarize in 
Section 6. 

2. DATASET 

In this section we describe the observations used to 
derive the extinction curve towards the Galactic Cen- 
ter. We used a VLA radio (2cm continuum) map, ISO- 
SWS spectra (2.4 to 45 ^m) and SINFONI imaging spec- 
troscopy (1.2 to 2.4 /im), of which the latter is used to 
construct line maps at Brackett-7 (2.166 /im), Brackctt- 
C (1.736 /im), and Paschen-/? (1.283 fim). The fields of 
view for the different datasets are compared in Figure [H 

2.1. Radio 

Figured] presents the radio and near-infrared data used 
combined with the field of views of all measurements. 
The radio data consists of combined 2 cm continuum ob- 
servations in the A, B, C and D configurations of the 
VLA. The phase center of the multi-configuration data 
set at 2cm is on Sgr A*. Standard calibration was done 
on each data set before each data set was combined in 
the uv plane. Details of phase and amplitude calibra- 
tion of each data set taken in different configurations of 
the VLA can be found in ,Yuscf-Zadeh fc Wardlc (,1993t ) 



Fig. 2. — Field of view for different measurements. In red, the 
radio 2 cm continuum image; in blue and green, respectively, are 
the H and Ks broadband images constructed from SINFONI HK 
data. The radio image is larger than displayed here, but the outer 
areas are not used in this publication. Light green lines mark 
the borders of the SINFONI field which has a hole around IRS7. 
Blue lines indicate the contours of the ISO PSF around 2.6 ^m, 
representing 90% , 50% and 10% of the central intensity. The red 
lines indicate the same for the ISO PSF around 19 /im. 



and' Yusef-Zadeh et al.l ()1998l ). The combination of these 
configurations makes high resolution possible, while at 
the same time ensuring good coverage of extended struc- 
ture. Due to the use of the smallest, D, configuration, 
we are able to detect structures of up to 50" in size. 
The resolution of the final image is 0.42" x 0.3" with the 
longer axis of the Gaussian oriented approximately north 
(P.A.« 4°). The image extends at least 48" from Sgr A* 
in all directions. We use a radio continuum map because, 
for the data available to us, the SNR and the resolution 
are better than for line maps. 

We chose to use a map at 2 cm because non- 
thermal emission, redder than the thermal emission of 
the minispiral, is les s dominant at shorter wavelengths 
(jScoville et al.l (|2003l ). on the other hand, used a 6 cm 
map) . Naturally, some non-thermal contamination is ex- 
pected at 2 cm as well. We control for this by testing 
whether there is a correlation between minispiral emis- 
sion and extinction strength (see Section 14. 2p . Sgr A* 
is a clearly visible non-thermal source in the radio map, 
see Figure [21 However, due to our high resolution it is 
well isolated from the minispiral. We are thus able to fit 
the source by a Gaussian and subtract it, smoothing the 
area around it. As a result, Sgr A* is not visible in our 
final radio map (Figure [3]). 

Two particula rly problemati c regio ns a re IRS2 and 
IRS13, for which [Roberts et al.l (|1996l ) and lShukla et al.l 
(|2004f ) derived significantly higher than average electron 
temperatures. In our analysis we simply ignore these 
regions, masking out IRS13 and IRS2. 

Difficulties can also arise, for interferometric observa- 
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Fig. 3. — Top: 2 cm continuum radio image from VLA. We 
subtracted Sgr A* by fitting a Gaussian to it. Tfie resolution is 
smoothed to 0.85". Bottom: Bradsett-7 map derived from SIN- 
FONI. Areas with stellar lines are replaced by the neighboring gas 
emission. Some areas at the border were not mapped and arc given 
zero flux. We fill the area around IRS7 that was not covered by the 
observations with the median flux of the surrounding. Both maps 
show the same area, with 250 mas pixel" ^, the same resolution, 
and the same color scaling. 

tions, from the fact that only a finite range of spatial fre- 
quencies can be sampled. This has the consequence that 
structure larger than the fringe spacing of the shortest 
baseline is not measured, which can result in an underes- 
timate of the flux, if the observed object is much larger 
than the largest angular scale sampled. For us, however, 
this is not a problem, since our map, due to the use of 
D configuration as mentioned above, we are able to de- 
tect structures of sizes up to 50" , while nearly all of the 
flux of the minispiral is concentra ted within 50" diameter 
(Figure 1 of IScoville etUI (l200l l. 

2.2. ISO 

We use the ISO-SWS data of lLuti (fTOOl . These data 
are reduced with the SWS interactive analysis system 
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Fig. 4. — ISO-SWS spectrum of the GC. The hydrogen lines used 
for the extinction measurements are marked by colored vertical 
lines. 



(jWieprecht et al.l I1998D . The data consist of a single 
aperture spectrum extending from 2.4 to 45 /im, shown 
in Figure ID 

The field of view of SWS (Ide Graauw et al.l [T996al) 
is approximately rectangular; it can be seen in Fig- 
ure [2l We consider t he deviations from rectangularity 
(|Beintema et al.ll2003[ ) in SectionO In our observations 
the longer axis of the field of view is oriented north-south 
(P.A. « 0°). The field of view is approximately 14" x 20", 
shortward of 12 /im, and 14" x 27" between 12 and 28 /im. 
The observations are centered on Sgr A*. 

2.3. SINFONI 

For shorter wavelengths (1.2 to 2.4 /im) we use spec- 
tra obt ained with the integral f ield spectrometer S IN- 
FONI (jEisenhauer et al.l l2003at iBonnet et all l2003| ) at 
UT4 of the VLTEI. We use the following data: R+K- 
band (spectral resolution 1500), seeing limited cubes 
(FWHM=0.85") covering an area of about 37" x 30" 
around Sgr A* from 2003 April 9 with a spatial sampling 
of 250 mas pixel" ^ x 250 mas pixel^^. J-band (spectral 
resolution 2000), seeing limited cubes (FWHM=0.6"-1") 
of most areas of the minispiral obtained 2010 May 21 
and July 6, with a spatial sampling of 125 mas pixel"^ 
X 250 mas pix el" ^. We apply the standard data reduc- 
tion SPRED ( Abuter et al.l[200l: ISchreiber et al.l[200l 
for SINFONI data, including detector calibrations (such 
as bad pixel correction, flat-fielding, and distortion cor- 
rection) and cube reconstruction. The wavelength scale 
is calibrated with emission line lamps and finetuned with 
atmospheric OH lines. The HK data is atmosphere cor- 
rected and kept in ADU count units. No conversion in 
energy is applied. We flux calibrate the J-band data with 
standards which were observed directly after the science 
observations, as is described in more detail in the next 
section. The Brackett-7 map, one of three linemaps de- 
rived from the SINFONI observations (described in the 
next Section), can be seen in Figure |31 

5 We use SINFONI data from ESO programs 70.A-0029(A), 
I83.B-100(P) and I83.B-100(R). 
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3. CONSTRUCTION OF LINE MAPS AND FLUX 
CALIBRATION 

3.1. Line maps 

In this section we describe the construction of 
Brackett-7, Brackett-C, and Paschen-/? line maps from 
the SINFONI data, see Table [H In order to do this, 
we integrate over the channels which contain the neb- 
ular emission and subtract the average of the adjacent 
spectral channels on both sides as background. The line 
maps contain the following types of artifacts, unrelated 
to the minispiral emission: 

• Bad pixels, in the sense of a large deviation of a 
few pixels from their neighbors. We identify them 
mostly manually. 

• Gaps between cubes, when not well aligned 

• Emission and absorption line stars 

We interpolate these artifacts, treating them like bad 
pixels in the data reduction. We also replace the hole 
around IRS7 in our map with the median flux of the 
surrounding for obtaining line maps, see the Brackctt 7 
map in Figure [3l 

3.2. Flux calibration 

We flux calibrate the Brackett-7 map with a NACO 
Ks-image of the GC from April 29, 2006. On this day, a 
standard star was observed with a zero-point uncertainty 
of 0.06 mag. In order to translate this calibration to the 
SINFONI data, we extract a mock Ks-image (Fsinks) 
from the SINFONI cube. This is done by multiplying 
the atmosphere corrected SINFONI data slice by slice 
with factors which represent the product of atmospheric 
and NACO Ks-filter transmissior0 for every slice. The 
factors are scaled such that their integral is one. In this 
way we calibrate a single SINFONI slice at the isophotal 
wavelength of the Ks-band. 

We measure the full flux in the NACO image and in 
FsiN Ks over the area covered by the SINFONI data. We 
estimate the uncertainty of the cross-calibration by di- 
viding the two images in nine parts and measuring in 
each part the count ratio. We obtain an uncertainty of 
0.03 mag from the standard deviation of the nine count 
ratios. 

In total we use the following factor for calibration of 
the Brackett 7 line map: 

j_ ^AKsAisoKs -|^Q-0.4ZP j-^-j 
i? ABr7 

^ Fx Ks Also Ks (zP.Aco+2.5 log t^^) 

i? ABr7 

Here, Fa ks is the Ks calibration of Vega 
(jTokunaga fc Vaccal l2005n . R the number of shces 
per iim and ZP the zero point in magnitudes of the 
SINFONI data. We assume an error of 0.03 mag 
for F\Ks- Because we use ADU spectral data, the 
amount of energy per ADU depends on wavelength. 

® We use the NACO Ks-filter transmission from the NACO web 
site: http: / / ww w.eso.org/ sci /facilities / par anal / instruments / naco | 
/ inst/tilters.html 




Fig. 5. — Color image from J-band SINFONI data: blue 1.18 /im 
image, green 1.3 fim image, red Paschen-/? line map. 



Since the calibration is a t the isophotal wavelength 
(|Tokunaga fc Vaccal I2005D . we multiply F\ks by 
AisoKs/ABr7- Compared to using a standard star, our 
calibration has the advantage that we use the same 
NACO data as ISchodel et "aIT(|2010D for calibration. In 
this way we reduce the uncertainty when comparing our 
results with theirs, see Section [521 

We calibrate Brackett-C in a similar way to Brackett- 
7, using the integrated flux of an H-band NACO image 
from April 29, 2006, which is compared to a H-band im- 
age constructed from the SINFONI cube. The errors of 
the calibration are 0.05 mag (NACO zero point), 0.04 
mag (cross-calibration) and 0.03 mag uncertainty from 
the conversion to energy {F\h)- 

We use data cubes flux calibrated with standard stars 
for Paschen-/? (see Figure O, because a flux calibrated 
GC image in the J-band was not readily available. By 
comparing the fluxes of different standard stars from the 
two nights in April and July 2010 in which the data was 
obtained, we validate that the calibration is stable, with 
only 0.034 mag rms scatter and 0.005 mag bias over the J- 
band spectrum. We test the calibration accuracy further 
by comparing adjacent areas in the data from April and 
July. In doing so, we do not find signs of significant dis- 
continuities in the Brackett-7 to Paschen-/? ratio at the 
border between the two areas. We estimate the system- 
atic calibration error to be the same as for Brackett-7: 
0.075 mag. 

For the ISO data we use the absolute calibration of 
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iLutz et alJ ()1996D . The calibration error dominates the 
total error for bright lines. It is 10% below 4.06 fim and 
20% above. 

Since the emission of the minispiral is not homoge- 
neous, it is necessary to compute the NIR and radio com- 
parison fluxes from the same area as the ISO-SWS beam 
profile. Because of this we mu ltiply the ISO-SWS beam 
profile (jBeintema et al.l I2003D in every ISO band with 
the Brackett-7 line map and obtain the Brackett-7 flux 
within the ISO-SWS beam profile. Since the Brackett-7 
map does not cover the full ISO beam, we also calculate 
the radio flux contained in the ISO-SWS beam profile 
in the same way. Over the Brackett-7 field the differ- 
ence in the beam correction is 3 % between using radio 
and using Brackett-7. This difference could be caused 
by non-thermal emission. Therefore, we use the radio 
data only for calculating how much fiux is missed due to 
the smaller field of the Brackett-7 data. In total, we find 
that the flux is about 7% smaller below 12 yum and about 
6% higher above 12 /im compared to a sharp 14" x 20" 
field. The sudden change at 12 /im is an artifact caused 
by the field size increase to 14" x 27" at this wavelength. 

4. DERIVING EXTINCTION FROM HYDROGEN LINES 

The ratios between different hydrogen recombination 
lines in HII regions depend only weakly on the local 
physics ()Hummer &: StorevI I1987D . Hence, it is possi- 
ble to derive the relative extinction (A) between two 
lines a and b by comparing the observed flux ratio 
(i^(a)obs/-F'(6)obs) with the expected (extinction-free) 
flux ratio {F{a)cxp/F{b)c^p): 



Aa-h 



-2.5 X log 



/ f(a)obs/F(6)obs 

U(«)cxp/i^(6) 



Jcxp 



(3) 



Using infrared data only, and without making further 
assumptions about the shape of the extinction curve, it 
is only possible to obtain the relative extinction. In order 
to obtain absolute extinction values it is necessary to use 
a wavelength for b at which the extinction is known inde- 
pendently, or negligible. Then it is possible to calculate 
the expected flux at a from the observed flux at b using 
the known flux ratios fo r emission from a gaseous nebula 
(jBaker k Menzellll938l ). We use their Case B, for which 
the nebula is opaque to Lyman radiation but transpar- 
ent to all other radiation (in Case A the nebula is also 
transparent to Lyman radiation.): 

i^(a)oxp = c X F(6)obs (4) 

c follows from the Case B calculation and depends on the 
radio frequency and infrared line used, as well as on the 
electron temperature, see also Section HTTl 

In our analysis, we use radio data for b and infrared 
lines in place of a. The formula for absolute extinction 
is thus: 



Air = -2.5 X log 



F{m)o 



bs 



(5) 



c X F(radio)obs 

We use the Case B line ratios in Hummer fc StorevI 
(|1987f ) for the Fobs ratios in Formula [31 The unextincted 
ratios of hydrogen emission originating from different 
ato mic le vels depend on the radiation state (like Case 
B) (|Lutz|[T999i) . As a first test of the vahdity of Case 



B we use the hydrogen line in the ISO-SWS spectrum 
at 7.50 fim, which is a blend of the lines 6-8 and 8-11. 
After accounting for the blending the 7.50 /im extinction 
differs only by 0.05 mag from the extinction for the 6-5 
line at 7.46 /im, assuming Case B. This implies that Case 
B is indeed valid for the Galactic Center ( Lutz 19991 ). 

4.1. Electron Temperature 

The physical conditions of the plasma, in particular 
the electron temperature (Te), have an important in- 
fluence on the extinction free flux ratios of the IR line 
fluxes to the radio flux . For example the dependence of 
Paschen a to radio continuum (free-free) emission (Sg) 
is Fpaa/S'ff oc (jScoville et al.l[2b03l Accordingly, 

it is necessary to know Tg for deriving the absolute ex- 
tinction. 

The electro n temperature h a s been derived at H92q! 
(8.3 GHz) by 'Roberts & Goss| (fl993[) and at H41a (92 
GHz) by Shukla ct al. (200^^ Both obtain T^ « 7000 
K in most parts of the minispiral. The consistency of 
both me asurements shows t hat all conditions, thus also 
Case B, ([Roberts et al.lll99lD for deriving Tg are fulfilled, 
even at the smaller and more problematic frequency of 
8.3 GHz. In this work, we take the measurement by 
[Roberts & Goss (1993) of Te = 7000 ± 500 K, since it 
has the highest SNR, however we correct it to account 
for the He"^ fra ction (see below) . 

According to [Roberts fc Gosa ([1993[ ). there is no sig- 
nificant spatial variation in the electron temperature. 
Thus, we think it is justified to simply assume a constant 
electron temperature. Even if the electron temperature 
varies, the variation does not matter as long as the value 
used is equal to the flux weighted average of Te of the 
area over which we integrate. 

The absolute value o f Te also depends on the He+ frac- 
tion: Y+ = He+/H+. [Roberts fc Gosi ([1993) derived a 
2cr limit for Y + of 3%, and measur ed Y+ = 5% ± 2% in 
another area. [Krabbe etlH ([199 If ) derived y+ « 4 % 
from the He I line at 2.06 /xm. We obtain Y+ = 2% ±0.7 
% using the same line. Considering all information, we 
assume — 3% ± 1 % for the calculation of Te. 

The elec tron temperature depe nds, according to 
[Roberts eF al. (1996) and .Roelfsema et al.[ ([1992D . in the 
following way on Y+: 



Te - T, 



1 



l + y^ 



0.87 



(6) 



Teraw IS the Te obt ained assuming = 0. Since 
[Roberts fc Gosd ([1993i ) used Y+ = 0, they calculated 
T e raw We therefor e calcu latc the real Te from the Te raw 
of [Roberts fc"Gosi (fl993l) using our value of Y~^ = 3%. 
We obtain Te = 6800 ± 500 K, and use this value here- 
after. The error in the electron temperature results in an 
absolute extinction error of 0.043 mag. In principle, it is 
also necessary to know the electron density of the emit- 
ting plasma to calculate the a bsolute extinction from in - 
frared lines and radio emission ([Hummer fc Storevll987l). 
We u se an electron density of lO^cm"'^ ([Shukla et al.l 
[20041 ) ■ but the results are not very sens itive to the ac- 
tual value used ([Hummer fc Storevl[l987l) . 

4.2. Extinction calculation 
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We first compute the Brackett-7 extinction from di- 
rect comparison of the Brackett-7 map with the radio 
data (Figure |3]), following Formula [5] We then calculate 
all other extinction values with respect to the Brackett-7 
data. This method has the advantage that we can se- 
lect the fiducial area by comparing with the high SNR 
Brackett-7 data (Table [T|). This is done by calculating 
the extinction relative to Brackett-7 extinction following 
Formula 131 The Brackett-7 extinction is then added to 
this value to obtain the absolute extinction. 

To calc ulate the expe c ted fl ux ratio we convert For- 
mula 3 of lScoville et all (|2003f ) to Brackett -7, our con- 
tinuum frequen cy of 15 GHz, and mJy, us ing the relative 
line strength in lHummer fc StorevI ()1987f ): 



[^P^] = (1.327 ±0.052) 10""ergs"^ cm"^ ^jy-i 

V ^ff J oxp 

Inserting this expected ratio and the measured 
Brackett-7 to radio ratio in Formula [5] we obtain a 
Brackett-7 extinction map. 

In order to estimate the impact of non-thermal emis- 
sion, we measure the resulting extinction in Brackett-7 
flux bins, see Figure El For most line fluxes, the ex- 
tinction is approximately constant. However for line 
fluxes smaller than 3.7 x 10~^^ ergss"^ cm~^ pixeP^ the 
extinction is anticorrelated with the Brackett-7 flux. 
Most likely the anticorrelation is caused by non-thermal 
emission of Sgr A East which becomes more impor- 
tant in fainter regions of the minispiral. Because of 
this, we mask out areas with fluxes smaller than 3.7 x 
10"""^^ ergss""'^ cm~^ pixel"^. The choice of this flux cut 
is somewhat arbitrary. However, since we only use the 
integrated flux of the whole ISO beam, which is domi- 
nated by the brighter fluxes, the uncertainty of the flux 
cut only introduces an extinction error of 0.03 mag. In 
addition, we mask out IRS2 and IRS13 (see Section [2TT|) . 
as well as IRS16NE and IRS16C because of remaining 
stellar emission line flux there. 

The extinction in the areas defined by t he above con- 
strain ts is similar to the extinction map of iSchodel et al.l 
()2010() . However, there is a scatter of about 0.4 mag in 
our map, when comparing pixel by pixel, due to the low 
SNR per pixel in Brackett-7 and radio data. Because of 
this, our extinction map is not of direct use. Neverthe- 
less, the small SNR per pixel does not affect our abso- 
lute measurement, because we use the flux from a total 
of 2854 pixels for derivi ng the Brackett-7 ext inction. 

For a comparison with'Schodel et ajj ()2010D we add up 
the Brackett-7 and radio fluxes wit hin 0.1 mag b ins. We 
use the extinction map of ISchodeleTaD (pOlOh to bin 
according to extinction. We then calculate the extinc- 
tion in every bin from the total fluxes therein by means 
of Formula [71 (see Figure 171). The differences between 
the extinction values from ISchodel et al.l ()2010[ ) and this 
work have an rms scatter of 0.10 mag. This scatter is 
bigger than would be expected from the SNR of both 
maps, but sy stematic prob l ems, s uch as residual aniso- 
planatism in ISchodel et al.l (|2010[ ). or Te variations for 
our map, could cause increased scatter. We assume that 
half of the scatter is caused by our data and add therefore 
0.07 mag as additional extinction error. This estimate is 




Brackett 7 flux [ergs/(s*cm^*pixel)] 

Fig. 6. — Extinction in the central 14" X 20", binned in Brackett- 
7 flux. The squares and error bars mark the median extinction 
per flux bin. (The error is based on the median deviation, scaled 
to a Gaussian Icr error by multiplying it with 1.483.) While for 
small fluxes the extinction is anticorrelated with the flux, there 
is no such correlation at fluxes larger than the vertical line at 
3.7 X 10~^® ergs s~^ cm~^ pixel"^ . The red line is the integrated 
extinction using all pixels with a flux larger than or equal than the 
flux in the given bin. There is nearly no variation. 
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Fig. 7. — Comparison of the exti nction derived from our data 
and that of the Ks extinction m ap of[Schodel et al.l (12 0101 ) . We bin 
our data according to the map ofLSchodcl ct al. (2010) and measure 
the extinction in every bin. The error bars show the expected error 
due to the single pixel SNR. The rms scatter in the plot is 0.1 mag, 
exceeding the expectation. 

conservative, because many problems should average out 
over the full area. 

Our aim is to measure the wavelength dependent ex- 
tinction curve in one line of sight. Spatially variable ex- 
tinction can however cause deviation of the line of sight 
extinction curve, if different measurement areas are used 
at different wavelengths. In addition, for comparison 
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TABLE 1 

Combinations of lines or continuum used for deriving extinctions 



Line A Lino B area used ('good' area) 



2 cm eontinuum Brackett-7 'ISO beam', excluding: low Br-7 flux, IRS16C, IRS16NE, IRS13, IRS2 

Brackett-7 Brackett-C 'ISO beam', excluding: low Br-7 Aux, Br-7/ Br-C outliers, IRS16C, IRS16NE, IRS13, IRS2 

Brackett-7 Paschen-/3 'ISO beam' inside the J-band data, excluding: low Br-7 Aux, IRS16C, IRS16NE, IRS13, IRS2 

Brackett-7 ISO lines in each case the adequate ISO beam 



with iSchodel et afl (|201Cll ). it is also necessary to correct 
for the difference between the average extinction of this 
map and the Brackett-7 flux weighted extinction in the 
same area. In order to take both effects into account, we 
calculate the averag e extinction of the extinction map of 
ISchodel et al.l (|2010[ ) within the region of the ISO beam, 
and the flux weighted extinction within the ISO beam 
for the different line maps. The difference in extinction 
is small for all lines. We find the biggest deviation is 
-0.03 mag for Paschen-/?, for which line the smallest area 
was covered, see Figure [5j We apply these corrections to 
our raw extinction values. 

There is an additional bias in integrated extinction 
measurements that can cause apparent flattening of the 
extinction law, and is especially irn portant for integrated 
galaxy SEDs ( Cal zetti et al.ll2000[ ). This bias arises be- 
cause the measurement of the extinction is flux-weighted: 
this flux being the observed (extincted) flux. This creates 
a problem, because if the extinction is spatially variable, 
the extincted flux will be brighter in low-extinction sub- 
regions of the image. Thus more weight will be given 
to these bright, low-extinction sub-regions of the image. 
A similar effect of very inhomogeneously distributed ex- 
tinction in the case of galaxy integrated measurements 
( mixed cas e ) is v isible in the top left panel of Figure 2 
in iCalzettil (|2001l ) . The effect of this bias is negligible 
for negligible extinction (i.e. at long wavelengths), but 
increases towards shorter wavelength together with the 
extinction, producing an apparent flattening of the mea- 
sured extinction law. We test and correct for this bias 
in our measurements, which is however small (up to 0.08 
mag, see Appendix |A| . 

In the H-band, the hydrogen lines are weaker, both 
intrinsically as well as being subject to higher extinc- 
tion. The atmospheric OH-lines are also stronger in H- 
band. Because of this, only Brackett-^ at 1.736 /im is 
easily detectable, although even this line is polluted by 
an atmospheric OH-line, despite sky subtraction. The 
structure of this remaining OH-line emission/absorption 
follows the cubes which are combined to the mosaic and 
is caused by the time variability of the sky. The sky 
has to be observed offset from the GC in space and thus 
also in time due to the brightness of the GC. To cor- 
rect the sky remnants, we derive OH-line strength maps 
from the next two strongest OH-lines at longer wave- 
lengths, and subtract the average of the two maps from 
the raw Brackett-^ map. We scale the subtraction such 
that any observable structure due to individual cubes 
in the combined cube vanishes. We use the uncertainty 
of the scaling factor and the difference between the two 
OH-line maps to estimate the error. In order to further 
exclude unphysical outliers we do not use pixels which 
deviate by more than 2.5 cr from the median Brackett-7/ 



Brackctt-C flux ratio. Integrating over the useful area of 
both lines we obtain FBr^y/FBrC = 16.1 ± 1.3. The error 
due to the OH line subtraction and the pixel selection 
uncertainty adds an error of 0.081 mag to the extinction 
for Brackett-C. 

For Paschen- /3 we used the smaller area (Figure [5]) 
within its overlap with the central 14" x20". We use the 
hydrogen lines detected in the ISO-SWS spectrum (Fig- 
ure |4]) to derive extinction values at the corresponding 
wavelengths. For the ISO data, which is not spatially 
resolved, stellar emission could be a problem. We test 
this by comparing the total Brackett 7 from our map 
with that computed upon declaring stars visible in the 
linemap to be "bad pixels" . The flux difference between 
the two cases is only 0.04 mag. This contribution, com- 
pared to the extinction error of the ISO data of at least 
0.14 mag, is therefore not relevant. Since the intrinsic 
(non extincted) line ratios follow Case B (Section 14. 2p . 
the intrinsic ratios are well known and do not cause an 
additional error apart of the Te uncertainty. 

4.3. Extinction derived 

The extinction derived in this analysis is presented in 
Table [2] and Figure [8] It decreases with wavelength, fol- 
lowing a linear relation in this log-log plot between 1.2 
and 2.8 /im. At longer wavelengths, the extinction curve 
changes shape: it contains more bumpy features, and 
is higher than expected from linear extrapolation. For 
the purpose of investigating the extinction, it is useful 
to separate the two regimes. Therefore, we define the 
near infrared (NIR) as the wavelength regime from 1.2 
to 2.8 fim and the mid infrared (MIR) as the wavelength 
regime from 2.8 to 26 fim. The linear relation within the 
NIR in the log-log plot implies a power-law: 

At^AaX {Xa/XbT (8) 

For the five NIR lines we obtain a = —2.11 ± 0.06 and 
ABr-y = 2.62 ±0.11, see Tabled The of the fit is 2.70 
given 3 degrees of freedom. Hence, a power-law is a good 
description of the NIR extinction. Since the extinction 
is obviously grayer beyond 3.7 //m, we test if there is 
any indication for fiattening within the NIR by fitting 
the red and blue part of the 5 NIR lines separately, see 
TableEl The slope of the red part, with a = -1.76±0.39, 
is only 0.9 a flatter than the blue part. Because this 
change in slope is not significant, we use the power- law 
obtained from all 5 NIR lines for deriving the broadband 
extinction values (Appendix [C]). The extinction between 
3.7 and 8 /im is fitted by the following power- law: 

A{X) = (1.01 ±0.08) X (A/4.9Aim)-°-4^±°-29 

Thus, the MIR extinction is grayer than the NIR extinc- 
tion. 
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TABLE 2 

Average hydrogen line extinction values towards 

THE central 14" X 20" OF THE GC 



A [fim] 


extinction 


1.282 


7.91 ± 0.11 


1.736 


4.30 ± 0.13 


2.166 


2.49 ± 0.11 


2.625 


1.83 ± 0.13 


2.758 


1.51 ± 0.19 


2.873 


1.84 ± 0.19 


3.039 


2.07 ± 0.19 


3.297 


1.66 ± 0.15 


3.74 


1.19 ± 0.14 


3.819 


1.19 ± 0.31 


3.907 


1.09 ± 0.3 


4.052 


1.01 ± 0.13 


4.376 


1.09 ± 0.29 


5.128 


0.99 ± 0.24 


5.908 


1.04 ± 0.24 


6.772 


0.84 ± 0.29 


7.459 


0.81 ± 0.23 


7.502 


0.79 ± 0.24 


8.76 


2.04 ± 0.34 


12.371 


1.34 ± 0.32 


19.062 


1.34 ± 0.5 





TABLE 3 




NIR infrared power-law extinction fit parameters 


Used range [fJ-m] 


Asr-f slope a 


xVd.o.l 


1.282 to 2.166 


2.60 ±0.11 -2.13 ±0.08 


2.05/1 


2.166 to 2.758 


2.50 ±0.11 -1.76 ±0.39 


0.54/1 


1.282 to 2.758 


2.62 ±0.11 -2.11 ±0.06 


2.70/3 



5. DISCUSSION 
5.1. Comparison with literature for the GC 

We now compare and discuss the extinction curve ob- 
tained in this work with the h terature. 

First of all, we compare withlRiek e fc Lebofskvl ()1985l ). 
who used mainly data from GC stars to derive the in- 
frared extinction law up to 13 ^m. Their finding, how- 
ever, a single power-law of slope « —1.54 from J- to 
M-band, is in contrast with our results. The differences 
probably lie (i) in their assumption of an universal IR ex- 
tinction law, whereby Ry is determined only indirectly, 
and (ii) variability of the stars used. Firstly, for abso- 
lute calibration of the extinction law they used o Sco, 
which is also detectable in the optical , together with the 
relatio n E y k/Eb-v =2-744 from iSchultz fc Wiemeii 
(|1975f) and iSneden et al.l (|1978[ ). However, this relation 
was obtained outside of the GC, making their extinc- 
tion law vulnerable to line of sight variations, see Sec- 
tion l5.4l Secondly, they employ the Ey-M / Eb-v of the 
GC and the bolometric luminosity of IRS7 to estimate 
Ry = Av/Eb-v = 3.09. In the process they use mea- 
surements of three stars in the M-band, which is domi- 
nated by atmospheric emission, to estimate their lower 
limit on Ry. Not only are all of the GC stars used by 
iRieke fc Lebofskvl (|1985l ) supergiants, which are in gen- 
eral variable, bu t variability has, i n particular, been re- 
ported for IRS7 ()Blum et al.lll996[ ). Their indirectly de- 
termined Ry could introduce a systematic error in their 
extinction measurement. For example, if we reduce Ry 



to 2.98 in their calculation, we obtain from their data 
a = —2.04 between J and K. This is compatible with 
our a. 

In a later work, iRiekd ()1999( ) updated 
IRieke fc Leb"ofek^ (p8^ using NICMOS data (FllOM, 
F145M, F160W, F222M), again of stars in theCC region 
They used spectra of the IRS16 stars (jTamblvn et alj 
|1996[) for deriving absolute extinctions. The extinction 
curve they found is flatter in 'HK' (a = —1.58) than 
in 'JH' {a ~ —1.95), in contrast to our extinction law 
(Section I4.3P and to most other publications which 
derive a constant slope in the NIR (e.g, iDraind 1X989') . 
In order to test whether the NICMOS data generally 
contradict our NIR extinction law (Section 14.31) we 
us e stars which ha,ve N ICMOS magnitudes published 
in iMaillard et al.l ()2004[ ). We use only stars which 
have been spectro scopicall y ident i fied as early- type 
(jPaumard et al.l [2006: Bartko et al.l l2009[ ). Early- type 
stars have colors « in the NIR. We exclude IRS13E2 
and I RS13E4.0 from t he analysis because they are 
dusty (|Fritz et al.l I2O1O0 . Using the rescaled versions 
of the extinct i on m ap and the Ks-band magnitudes of 
iSchodel et all (|2010l ). the NICMOS data are consistent 
with a single power-law with a slope of a = —2. 16 ±0.08. 
Thus, the NICMOS data do not contradict our result of 
a = -2.11 ± 0.06 fS ectionlQl). 

In another study, iViehmann et al.l (|2005[ ) performed 
seeing-limited ISAAC L- and M-band photometry. As- 
suming stellar colors they derived a flat extinction slope 
from L to M: Am/Al = 0.966 ± 0.05. From our 
NACO broadband extinction values (Appendix [Cj we 
obtain a consistent value: Am/Al — 0.88 ± 0.23. How- 
ever, these resul t s are not directly comparable, because 
IViehmann et al.l (|2005l ) used a narrow M-band with a 
FWIIM= 0.10 /Ltm at 4.66 ^m. At this wavelength, 
a CO absorption feature is visible in the ISO spectra 
(iLutz et al.lll996HMoneti et al.l l200l1). see Figure p 

iScoville et al.l ()2003l ) obtained an absolute extinction 
map of the GC for Paschen-a. To do this, they used 
NICMOS narrowband imaging at Paschen-a (1.87 ^m) 
and at 1.90 /^m to construct the line map. They took 
a VLA continuum map at 6 cm to use as extinction- 
free data. Using their extinction map, we compute an 
average extinction of Apaa = 3.54 integrated over the 
ISO beam. This value is consistent with our value of 
Apaa = 3.56 ± 0.11, obtained from the best-flt power- 
law. 

ISchodel et afl ^M) measured the extinction towards 
the GC using NACO data in H-, Ks- and L'-band. By 
comparing the magnitude of the peak of the luminos- 
ity function (the red clump) with its expected magni- 
tude, they obtained the total light modulus of the GC. 
They then used a distance to the GC of 8.03 ± 0.15 
kpc, derived from different works, to obtain the extinc- 
tion: Kh = 4.48 ± 0.13, Kks = 2.54 ± 0.12 and 
A f/ = 1.27 ±0 . 18. 

ISchodel et al.l (|2010D also derived the selective extinc- 
tion Eh-Ks = {Ah — Aks)IAks from the red clump ex- 
tinction values. They then used Eh-Ks for obtaining a 
Ks-band extinction map from the observed stellar colors. 
Surprisingly, the average of the map is, with Aks = 2.70, 
larger than the extinction obtaine d from the red c lump . 
It can be seen from Table A.2 in ISchodel et all polol ) 
that most individual stars have larger values of H-Ks 
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TABLE 4 

Error sources for the extinction 



Error source 


Amag Amag 

for Brackett 7 for other lines 


• IR calibration error 






and line SNR 


0.073 


0.073 to 0.50 


• electron Temperature 


0.043 


0.043 


• selection of good pixels 






in Brackett-7 


0.03 


0.03 


• scatter comparing 






with Schodel et al. ('2010) 


0.07 


0.07 


• selection of good pixels 






in other lines 





to 0.088 


total 


0.114 


0.114 to 0.507 



than the stellar population of the red clump in the lumi- 
nosity function, the reason for whi ch is not clear. 

We do not know which of the two lSchodel et al.l (|201Clf ) 
extinction values should be preferred. Because of this, we 
use the average and enlarge the error by adding ^1/2 x 
iAKsi-AKs2) = 0.11 mag. We obtain = 2.62±0.16 
for the full field of lSchodel et al.l (^2,010). Their extinction 
map has an average of Aks = 2.70 over the full map. 
However, the average is Aks = 2.68 over the ISO field of 
view. Therefore, we correct the average by multiplying 
with 2.68/2.70 and obtain Aks = 2.60 ± 0.16 as a final 
value for comparison with our extinction. 

The issue of two possible extinction values in the Ks- 
band also extends to the other bands. We assume a 
linear scaling of the extinction values based on stellar 
colors with the extinction values derived from the red 
clump in each band. We again use the average of red 
clump extinction and the selective extinction. We ob- 
tain = 4.58 ± 0.24 and = 1.30 ± 0.19 as final 
values for comparison with our extinction values. 

If we now compare with our broadband extinction 
values (see Appendix [UJ, we obt ain the follow i ng dif - 
ferences between our results and ISchodel et al.l (|2010[ ): 

linos — AffSch = 0.07 ± 0.25, lines " ^if s Sell = 

0.07±0.17and^L/ii„es-^L'Seh = -0.10±0.23. For the 
calculation of the errors here, we exclude the zero-point 
errors in the H and Ks-band because both results use the 
same calibration data, and we also exclude the error due 
to Ro because we use the differences to measure Rq , see 
Section 15.21 To su mmarize, our values are all consistent 
with the values of iSchodel et al.l ()2010[) . 

In Table [J] we present the various sources of error 
for our extinction values. Our data are well fitted 
with a power-law in the NIR {x^/d.o.f — 0.71), see 
Section 14.3! as expected from the established literature 
(fCard elli et al. 1989.). An error in the radio data (such 
as an error on Tg) affects all values in the same way, 
such that its influence on the of power-law fit is 
much smaller than the infiuence of independent errors 
in the NIR data. Because the power-law does not con- 
tinue to the MIR till to negligible extinction in our data, 
it is not possible to obtain the absolute e xtinction inde- 
pende nt of the radio data, as was done in iLandini et al.l 
(| 19841 ) . Howev er, the value for T^. whic h we use is consis- 
tent w ith both lRoberts fc GossI (|1993( ) and lShukla et al.l 
(j2004l ). Therefore, it is unlikely that there is a relevant 
systematic error in the extinction due to the radio data. 

5.2. A photometric distance to the GC 



The distance to the Galactic Center Rq is used to 
derive distances to all but the closest regions of the 
Galaxy and is one of the fundamental parameters for 
building Milky Way models. Knowing Rq is hence of 
general relevance. Prior to this publication there was 
no reliable measurement for the GC extinction that was 
independen t of Rq. Works that d i d derive Rq photo- 
metr i cally (iNishivama et all l2006bt iGroenewegen et alj 
l2008t iDambisI 120091: iMatsunaga et al. I I2009D used stars 
in the bulge, not directly in the GC. On the other 
hand, stars directly in the GC have been used to de- 
rive Ro from the dynamics of their st e llar orbits around 
the SMBH (lEisenhauereFaLl [2003bl: iGhez et al.l [20081: 
iGillessen et al.l I2009D , and by using the statistical par- 
allax of th e populatio n of late-type stars in the GC 
()Genzel et al. 2000; Tri ppe et al ■1120081 ) . 

With our extinction measurement we can obtain a pho- 
tometric Ro from stars in the GC. We combine our ex- 
ti nction measurement with the extinction measurement 
of lSchodel et al.l (|2010[ ). These authors obtained the to- 
tal luminosity modulus of the red clump stars in the GC. 
Because the luminosity modulus involves both the ex- 
tinction and the dis tance modulus it was necessary for 
ISchodel et al.l ()2010[ ) to assume a distance for calculating 
extinction values, for which they used Rq = 8.03 ± 0.15 
kpc. Here, we can use the extinction differences between 
the two works, see Section l5.ll in order to estimate Rq. 
For the Ks-band, we obtain Rq = 7.78 ± 0.63 kpc. From 
the H-band we obtain Rq = 7.78 ± 0.95 kpc, and from 
the L'-band we obtain Rq = 8.41 ± 0.94 kpc. The er- 
rors follow from the errors of the extincti on differences 
betwe en our extinction values and those of lSchodel et al.l 
()2010l) , see Section 15.11 We use the weighted average of 
all values as our final value: Rq = 7.94 ± 0.65 kpc. For 
the error we use the smallest single error, the relative er- 
ror of Ks-band, because the errors are correlated between 
t he filters. 

INishivama et al.l (|2006bl ) used red clump stars in the 
inner bulge to derive a photometric distance and ob- 
tained Ro = 7.52 ± 0.36 kpc. Our value is about 0.6 a 
larger. Since in both works the same absolute red clump 
magnitude is used, the consistency indicates only that 
other uncertainties (like the extinction) are not larger 
than assumed. However, it is still possible that the 
absolute magnitude has a bigger error than ass umed. 
Our result is also consistent with the review of iReidI 
11993), who det ermined Rq — 8.0 ± 0.5 kpc, and the 
recent review of iGenzel et al.l (|2010l ) , who determined 
Ro = 8.15 ± 0.14 ± 0.35 kpc from direct and indirect 
measurements, and Ro = 8.23 ± 0.2 ± 0.19 kpc from only 
direct measurements. 

Conversely , assum ing that the direct estimate of Rq by 
IGenzel et all (l2010l ) is correct we can test the red clump 
magnitude: The red clump in the GC has then a magni- 
tude of Mks = —1.5 9 ± 0.13. The value used by us and 
ISchodel et all (|2010[ ) is Mks = -1.47, and the two values 
are consistent, alth ough the GC star formation history 
(jBlum et al. (2003), Pfuhl et al. in prep) and metallic- 
ity (Cunha ct al. 2001) was not model ed to obtain Mks- 
There fore, as expected according to iSalaris fc Girardil 
(|2002D . Mkb of the red clump is therefore relatively in- 
dependent of the star formation history and metallicity 
and thus is a reliable distance indicator. 
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5.3. Spatial distribution of the extinction and the 
extinction towards Sgr A * 

Our data do not have enough SNR to obtain a good 
extinction map. However, because our absolute values 
fo r the extinction are m ore accurate than the extinction 
of iSchodel et all (|2010D , combining our abso lute value 
and the extinction map of ISchodel et al.l (1201 Of) is useful. 
For this, we adjust the extinction map of ISchodel et al.l 
([2010) such that in it the extinction is the same as our 
interpolated Asr 7 of t he same area . This means that we 
multiply the map of iSchodel et al.l (|2010[ ) by 0.976. 

From the adjusted map we use Ksr-y towards Sgr A* 
to derive broadband extinctions (Appendix [C| of: Ah = 
4.21 ± 0.10, Aks = 2.42 ± 0.10 and Al' = 1.09 ± 
0.13. Effectively, this is the extinction to the stars with 
r< 0.5" around Sgr A* due to the procedure used by 
ISchodel et aII(|20li3D . However, it is adjusted for the SED 
difference between stars and Sgr A* (Appendix ICj) . 

The light that reaches Earth from the GC crosses many 
regions of the Galaxy. It is possible that dust associ- 
ated with the minispiral extincts the light in the GC 
itself. However, because the extinction derived from the 
minispir al is consistent with the extinction derived from 
stars of ISchodel et al.l ()2010f ) , the extinction must oc- 
cur mainly in front of the GC. For testing this further, 
we smooth the Br ackett-7 flu x to the resolution of the 
extinction map of Sch odel et al.l (| 2 pi0f). This we co m- 
pare with the extinction map of ISchodel et al.l (|2010l ) in 
bins defined by the smoothed minispiral flux. We do not 
find a correlation between the minispiral flux and the ex- 
tinction derived from stars. Furthermore, the scatter of 
the median extinction is only 0.04 mag over the differ- 
ent flux means. In addition, the small fa r infrared flux 
(jBecklin et al.lll982t iGuesten et allll987| ) of the central 
few parsecs shows that hardl y any UV radiation is ab- 
sorbed there (jBrown fc Lisztl 119841 ) . Ah in all, the NIR 
extinction must be very small (at least A^s < 0.1) inside 
the central parsec. 

In order to further constrain the location of the ex- 
tinction in the line of sight, we use the H — Ks values 
of the stars in Table A.2 of ISchodel eFal] (|2010D . Since 
the intrinsic color \H — Ks\ < 0.2 for nearly all stars, the 
H — Ks of each star depends nearly only on the extinc- 
tion. We then measure how many stars in the table have 
a H — Ks compatible with zero extinction. Thereby, we 
exclude stars which are so blue and faint in Ks that they 
would be too faint for detection if they would have the 
extinction of the GC. We exclude them in order to avoid 
a bias towards foreground stars. After this, we obtain 
that only five of 6324 stars have a H — Ks color com- 
patibl e with zero extinction. According to iPhilipp et al.l 
(|1999| ). the Galactic disk and bulge (outside of 300 pc) 
have 2.3 % of the flux of the GC at r= lO'Q Because the 
ratio of extinction free stars to all stars of 0.1% is much 
smaller than the fiux contribution of Galactic bulge and 
disk to the fiux in the center of 2.3 %, there must be ex- 
tinction within the Galactic disk. We then measure up 
to which extinction it is necessary to include stars in or- 
der to account for the 2.3 % star contribution of Galactic 

''' We exclude the central 10" because there the total flux is dom- 
inated by young stars. Since there are only a few young stars due 
to the top heayy IMF l|Bartko et al.|[2010l ) they are irrelevant for 
the star number ratio. 



TABLE 5 

NIR INFRARED EXTINCTION LAW FROM LITERATURE 



Publication 



i/Ei-, 



slope a 



Indebetouw et al. (2005') 
Messineo et al. (2005) 
Nishiyama et al. (2006a) 
Straizys & Laugalys (2008) 
Gosling et al. (2009) 
Nishiyama et al. (2009) 
Stead & Hoare (2009) 
Zasowski et al. (2009) 
Schodel et al. (2010) 
Our work 



2 ± 0.13 



2.09 ±0.13 



2.11 ± 0.1 



-1.65 ±0.12 
-1.9 ± 0.1 
-1.99 ±0.08 
-2.07 ±0.23 
-2.64 ±0.52 
-2.23 ±0.23 
-2.14 ±0.045 
-2.26 ±0.17 
-2.21 ±0.24 
-2.11 ±0.06 



weighted average 



-2.07 ±0.16 



Note. — NIR extinction slope measurements of diffuse ISM 
extinction since 2005. When a is not given in the pubUca- 
tion, or the effective wavelength arc unusual we calculate a from 



oft 



and 2.164 /im JNishivama et al.l2009l) . The error is either the mea- 
surement er ror in the publication o r the scatter of different sight 
lines. From iNishivama et al.l i200g|) we only use the 2MASS data , 
since the SIRIUS data are identical to INishivama et al.l ]2006al) . 

We use all values for the calculation of the weighted average. 

disk and bulge. We find that it is necessary to include 
stars with extinction up Aks = 2.0. This is about 3/4 of 
the total extinction towards the GC. 

Therefore, the measured extinction is mainly not as- 
sociated wit h the Giant Molec ular Clouds in the Nu- 
clear bulge (jMezger et al.l 119961 ) and as such is not re- 
lated to special processes in the central 50 pc. Because 
most bulges do not contain a lot of dust, it is likely that 
most of the extinction is caused by dust in the Galactic 
disk. As a result, the measured extinction curve is likely 
a typical extinction curve of dust in the Galactic disk. 

The extinction towards the GC is higher than the 
average extinc tion of the bulge behind the Galactic 
plane (M arshall et al.l [2006) . However there are also 
regions with much h igher extinction close to the GC 
(|Ramirez et al.ll2008[ ). Therefore, the extinction of the 
Galactic Center is not exceptional for an 8 kpc view 
through the Galactic disk. 

5.4. The NIR extinction 

Our NIR data can be well fitted with a power-law of 
a = -2.11 ± 0.06 {x^/d.o.f. = 2.70/3). This strength- 
ens the case for the use of a po wer law as model for 
the extinction in the NIR (e.g., ICardelli et a l. 1989). 
Our slope is steeper t han the slope of a k, —1.75 of 
most revie ws, see e . g^ iSavage fc MathisI (|1979[ ): iMathii 
(fT990l) and lDraind (I2003D. However, most of the mea- 
surements (e.g.. 'Schultz fc Wiemeii Il975t iLandini et al.l 
[i984; Whittet 1988; He et al. 199^combined in these re- 
views used relatively few stars; these stars are detectable 
in the optical and partly also in the UV, have Ay < 5 an d 
are mostly closer than 3 kpc, see e.g. iHe et al.r(ll995D . 
The extinction measurement of lRieke fc Lebofskvl ()1985| ) 
towards the GC was also tied partly to measurements of 
stars which are visible also in the optical. 

Since about 2005, large infrared surveys have become 
available and are now used by most publications about 
extinction. Since now no detection in the optical is nec- 
essary, and because it is easier to characterize high ex- 
tinction, many publications (see Table [5]) measured the 
extinction towards the highly extincted inner Galactic 
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disk and bulge, at about 8 kpc distance. Most of these 
publications measure a ~ —2.1. 

It is possible that the change in the me asured a aroun d 
2005 is due to systematic errors. Stead fc Hoard (120091 ). 
for example, suggested that the use of the isophotal wave- 
length instead of the effective wavelength caused the flat- 
ter slope in the measurement prior to 2005. We think, 
however, that it is unlikely that errors in the effective 
wavelengths are the reason for the discrepancy: while it 
is correct that using the isophotal filter wavelengths can 
lead to errors jn a,_using the effective wavelength as pre- 
sented in lStead fc Hoard (|2009( ) overestimates a slightly, 
see Appendix |D] Furthermore, even if the isophotal 
wavelength is used, the systematic error on a is at maxi- 
mum 0.07 for the hot stars (see Appendix [Dl for a 9480 K 
star) used i n mo st works before 2005. Additionally, 
iFitzpatriclj ()2004l ) used 2MASS data of solar neighbor- 
hood stars and, using the effective wavelength, obtained 
a « -1.84. 

Since the dust probed by studies which obtain a ~ 
— 1.75 and a « —2.1, respectively, is not identical 
it seems likely that the extinction law varies between 
these regions. A strong piece of evidence for truly 
variable extinction i s the c orrelation of Ky with a in 
iFitzpatrick fc Massal ()2009D . whereas for the standard 
Ry = 3.1 a = -1.77 ± 0.05. This work used 14 stars 
with observations from 120 nm to Ks-band. All in all, 
we think it is likely that there is a transition of a mostly 
flatter NIR extinction in the solar neighborhood to a 
steeper one in most parts of the Galactic disk. Natu- 
rally, further tests of the NIR extinction slope variation 
via measurement of the NIR extinction slope in the local 
low extinction sight lines, and in the Galactic disk, using 
the best methods available today would be very valuable. 
For most molecular clouds, the extinction is, with 
pa -1.8 (Roman-Zimiga et al. 2007; Naoi ct al. 2006; 



Kenvon et al.iil998;,Flahertv et al.,,2007;,Lombardi et al.i 



!nY( 

200a) . flatter than for the Galactic disk sight lines. How- 



ever, there are also clouds with deviating a, with an 
a range from -1 to -2.4 ( Froebrich fc del Burgol [20061 : 
IWhittetJ[T988tlRacca et al.ll2002l : iNaoi et al.ll2007n . 

There are some features typical of m olecular clouds 
that a re visible in the ISO-SWS spectrum. iWhittet et al.l 
(|1997f ) determined that about a third of the extinction 
towards the GC is caused by molecular clouds with sub- 
structure and the rest by diffuse extinction. A relatively 
small contribution of the extinction by molecular clouds 
is also supported by the fact, that the extinction vari- 
ation towards the central parsec is less than one third 
of the maximu m extinction, as can be seen in the ex- 
tinction map of iSchodel et al.l (|20100 . This means also 
that the molecular cloud contribution is not visible in 
terms of star counts as it is in the case of e.g. the coal- 
sack. Accordingly, the line of sight towards the GC is 
not dominated by molecular cloud extinction. Thus, we 
exclude molecular cloud extinction in our quantitative 
comparison with the GC. 

For comparison we use publications about the extinc- 
tion towards the Galactic disk and bulge. In practice 
this means only publications since 2005, see Table El 
We use these data for calculating the weighted average: 
a = —2.07 lb 0.16. Of these 9 publications used, only 
llndebetouw et ahl ()2005f ) is inconsistent with the others. 
Hence, the extinction law is likely constant towards the 



inner Galactic disk and bulge and can probably be used 
also for other flelds in that region. 

Due to the high absolute extinction in the inner Galac- 
tic disk, this region contributes more to the global extinc- 
tion of the Galaxy tha n the solar nei ghborhood, or the 
halo of the Galaxy (.Schlaflv et al.l[20101 Therefore, the 
steeper extinction law in the NIR of the inner Galactic 
disk is probably more important when integrated over 
the volume of the Milky Way. Similarly, for other galax- 
ies for which our Galaxy is typical, a steeper extinction 
law could also be more important. 

5.5. The optical extinction 

iRieke et al.l (|1989[ ) estimated an optical extinction of 
Ky = 31 towards the GC using the extinction law of 
IRieke fc Lebofskvl ()1985f ). The GC is, due to this high ex- 
tinction, undetectable in the optical. As a result, any es- 
tim ate of Ay for the GC i s indirect. Since both this paper 
and ISchodel et al.l (I2010D obtain a diffe rent infrared ex- 
tinction law to IRieke fc Lebofek^ (fT985l) . it seems possi- 
ble that the optical extinction towards the GC is also dif - 
ferent from the one assumed in lRieke fc LebofskyI (|1985f ). 

In order to test this, we use data which has a more 
direct c onnection to the GC th an o Sco which was 
used by IRieke fc Lebofskvl (|1985f ). The shortest wave- 
length at which the GC has been observed is the z-b and 
(jHenrv et al.|[l98l lRosaeraI|[l992l: [Liu et al.lll993l) . In 
order to estimate the z-band extinction, we use the mag- 
nitudes of the IRS16 stars from these works, because 
these stars and their intrinsic c olors are well kii own. We 
neglect the data for IRS16C in lLiu et"aLl (|1993[ ) because 
the star is brighter in this work compared to all other 
works. W e use the extinction map and the IRS 16 mag- 
nitudes of ISchodel et al.l (|2010( ) in order to calculate the 
dereddened Ks-band magnitudes of the IRS16 stars. We 
then subtract the dereddened Ks-band magnitudes and 
an intrinsic color of z — Ks = —0.25 from the measured z- 
band magnitudes of IRS16 to obtain z-band extinctions. 

We convert the z-band extinctions into an extinc- 
tion power-law between Paschen-/? and the true effective 
wavelength for extinction measure ments (Atmp) of the 
GC (see Appendix [DII. In the case of lHenrv et all (jl984l ) 
and lLiu et al.l (119931 ) the wavelengths given in the publi- 
cations ar e probably A true ■ W e assign to them an error of 
0.01 /J,m. iRosa et al.l ([1991 give the central wavelength 
for stars extincted by GC extinction and unextincted 
stars. Since both are incorrect, see Appendix [Pl we cal- 
culate from the given central wavelength and FWHM of 
the optical system Atrue = 0.97 /zm. Due to the ad- 
ditional uncertainties in calculating the effective wave- 
length we assume an error of 0.015 /im. The effective 
wavelength uncertainty and the scatter of the IRS16 stars 
are used for the calculation of the errors. 

In this way we obtain a = 1.99 ± 0.09, a = 1.91 ± 
0.14 and g = 2 165 i b 0.13 from IHenrv et all ([1984); 
IRosa et al.l (|1992D and iLiu et all ([19931). The average 
slope is, with a = —2.02 ± 0.07, consistent with our de- 
termined NIR-slope of a = 2.11 ±0.06 f Section |43|). Us- 
ing the average of the three Paschen-/? z slopes we obtain 
Ai^m = 13.11 ± 0. 30. This ext i nction is slightly higher 
than interpolating IRieke et al.l ([1989D although our K- 
band extinction is smaller, see Figure [S] 

We now consider how the z-band extinction should be 
extrapolated into the optical. Molecular cloud features 
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are probably responsible for a third of the extinction to- 
wards the GC, see Section [5731 An above average con- 
tribution of molecular cloud extinction to the GC is also 
supported by the fact, that the extinction towards the 
GC is higher than the average bulge extinction (Sec- 
tion l5.3|) . If we assume a high Ry = 5.5 for the molecular 
cloud extinction and 3.1 for the ot her two thirds, t he av - 
erage Ry « 3.9 for the GC. Use a lCardelh et"al] (|l989f ) 
curve with this Ry (but with a = —1.85 between 0.91 
fiia and 1 fim for reducing the jumps in the slope) to 
extrapolate to the optical, W6 obtain Q 55^jn — 

30.3. 

However, most molecular clouds do not have an Ry 
of 5.5. Thus, even assuming that one third is caused 
by molecular cloud extinction and two thirds by nor- 
mal extinction, the real Ry is probably closer to 3.1. 
In addition, the high strength of aliphate (3.4 /im) and 
silicate (9.7 /im) features are different from the extinc- 
tion in molecular clouds (Section 15. 6p . These features 
are even stronger than in most diffuse extinction sight 
lines. Furthermore the steep NIR extinction slope of 
a = —2.11 is not typical of molecular cloud extinction 
(Section |5.4[) . Therefore, we obtain another estimate, ex- 
trapolating to the visible also with the normal Ry = 3.1 
(using a = —1.95 between 0.91 /im and 1 /im). In this 
case we obtain A_\— q 55^^! = 33. 

A much steeper extinction law than the s tandard Ry = 
3.1 fr o m B- b and to J-band is measu red f Udalski' '20031; 
Sumil 120041 : IRevnivtsev et al.l 120101: [Nishivama ct al. 
20081) towards parts of the bulge which have a much 
smaller extinction tha n the GC (Ay < 7 in case of 
INishiv ama et aTl ()2008f )'). In particular. iNishivama et al.l 
(120081) measured Ay /Aj = 5.32±0.14, while according 
to iRieke fc Lebofskvl (|l9 8 5) Ay/Aj = 3.55 ± 0.16 for 
the standard Ry = 3 . 1. W e fit these bulge data points 
with a CardcU i et al.l ()1989f ) curve with Ry — 2.0 (using 
a = —2.02 between 0.91 /im and 1 /im). Extrapolating 
the z-band extinction of the GC with this curve we ob- 
tain Aa^o.ss^hi = 44. The aliphate and silicate features 
towards the GC are stronger than in diffuse extinction, 
while they are even weaker in molecular cloud extinction, 
see Section 15.61 This means, that the strength of these 
features is anticorrelated with Ry for 3.1 < Ry < 5.5. 
Extrapolating this anticorrelation to the stronger feature 
towards the GC implies, that Ry < 3.1 towards the GC, 
which in tur n implies an Ay > 33. 

Figure 4 in Fitzpatrick & Massa] (|2009D shows a corre- 
lation between the infrared power- law slope a and Ry. 
The correlation implies that the slope in the infrared ap- 
proximately continues into the red. Using a second order 
polynomial to fit their 13 measurements (neglecting one 
outlier) we obtain Ry = 2.48 ± 0.06 for the GC, using 
our measurement of a — —2.11 ± 0.06 in the GC. Thus, 
because a small Ry implies a high Ay , this again implies 
the optical extinction towards the GC is probably large. 

X-rays can shed another light on Ay. X-ray 
photons, for example, are attenuated by scatter- 
ing and absorption which are related to the extinc- 
tion JMorrisonfc McCammonI [19831: IPredehl fc SchmittI 
[T99I . At the energy at which the GC is observed 
(« 3 to 8 keV), absorption by astronomical met- 
als is much strong er th an scattering by du st grains 
(jPredehl fc SchmittI [T995HPorauet et al.ll2008[ ). For the 
absorption the chemical state of the absorbing material 
is irrelevant and depends only on the integrated column 



density of astronomical metals. Thus the X-ray absorp- 
tion can be used to constrain dust models and the ex- 
tinction. 

Here we investigate the implications only for the ex- 
tinction. Observatio nally, the X-ray absorptio n is well 
correlated with Ay (jPredehl fc SchmittI 119951 ) in most 
Galactic sight lines. However, there a re AGNs with 
differ ent X-ray absorption to Ay ratios (jMaiolino et al.l 
120011: iLiI 120071) . Towards Sgr A* we obtain a column 
density of Nh = 10. 5 zt 1.4 x lO^^cm"^, averaging over 
the different states (jPorguet et al.ll2008l ). Thereby we 
assume an intrinsic X-ray powe r-law for the flares of 
Sgr A* ()Dodds-Eden et al.ll2009l ). We thus obtain an 
X-ray derived Ay = 56.7 ± 7.4 for the GC (Figure Et 
using the N///Ay relation of lPredehl fc Schmitt ('19951)'. 
Usi ng UV de rived /Ay relat ions (Bohlin ct al. 197^ 
iDraineHig^ iZubko et aLll2004D we obtain Ay from 53 
to 59. Our best fitting dust model (Section 15. 8p gives 
Ay w 48 for the measured Nh- 

There is the possibility, that part of the column den- 
sity towards Sgr A* a rises in ionized gas in the halo of 
the SNR Sgr A East (|Maeda et al.ll2002l: iPorguet et all 
|2008[ ). in which case not all of the Nh (measured by X- 
rays) towards Sgr A * is caused by dust. For example, 
iSakano et all (|2004D measured an Nh = 15 x lO^^cm"^ 
for hot plasma, and Nh = 7 x lO'^^cm^^ for somewhat 
colder plasma of Sgr A East. The second value is per- 
haps typical of plasma lying in the line of sight to Sgr A*. 
Accordingly, it is possible that the column towards the 
colder plasma corresponds to the column to Sgr A*. If 
this is th e case we obtain a lower Ay = 37.8, using the re- 
lation of P redehl fc SchmittI (|1995i ). However, probably 
the foreground extinction towards Sgr A East and Sgr A* 
is not fully homogeneous. In principle it is possible to ob- 
tain a better Nh towards the Sgr A region by measuring 
Nh towards point sources outside of Sgr A East. By 
comparing such a map with an IR excess map, it should 
be possible to obtain a better estimate for Nh for the 
GC. 

Overall, the evidence for an Ry < 3.1 is slightly larger 
than for Ry > 3.1, which means that, A;j— o.55/jm > 33. 
However, it is clear that only a direct measurement in the 
visible can clarify the value of Ry. The GC extinction 
attenuates 16NW, the brightest blue star in the GC with 
unextincted my « 5.5, to a magnitude of my ~ 34 to 
43 depending on Ry (since this simulated measurement 
integrates over the V-band the extinctions are smaller 
than the corresponding 0.55 fim extinctions). 

5.6. The MIR extinction 

At wavelengths > 2.8 /im extinction features are visible 
in the ISO-SWS spectrum, even shortward of the deep 
sihcate feat ure at 9.7 /im (Figure [4l). The se features ar e 
known fromlButchart et al. (1986); Willn er et al.l(|1979D : 
iLutz et al.l (|1996D and lChiar et al., (,20001 . 

Most of these featu res are caused by ices like the 
H2O feature at 3.1 tim (iButchart et al.l[T986l : iLutz et al.l 
[1996: Chiar et al.l l2000l) and the CO9 feature at 4.3 



/im (ILutz et all 1199 
othe rs, iWhittet et al 



de Graauw et al. 1996b2. Amoni 
11988,); .Rosenthal et all ([20a 



and iKnez et all (|2005D detected these features in many 
other sight lines with different strengths, both compared 
to one other as well as to the continuum extinction. 



IR extinction 



13 



o 

CO 



< 

o 
o 

.1— I 

!x! 
H 



CO 



• Porquet et al. 08 

• Bulge data 

X Sakano et al. 04 



•• 'Cardelli et al. 89' Ry=2.0 3.1 3.9 
- Cardelli et al. 89 Ry=3.1 



Rieke el al. 85/89 




♦ from continuum 
■ from hydrogen lines 



0.5 



10 



20 



X [/.im] 



Fig. 8. — Extinction towards the central 14" x20" of the GC. We use hydrogen lines for obtaining the extinction between 1.28 and 18 /^m 
(red boxes) and stellar colors for 1 /im (black diamond). We interpolate the data by use of the continuum emission (black line). The spectral 
resolutio n of the interpolation is not high enough to resolve all features like the CO feature at 4.7ftm fully. We mark the larger ex tinction 
featur es IILut z et al. 1996; Chiar et al. 2000). We use the central bulge data of Sumi (2004); Nishiyama ct al. (2008) and Revnivts cv et ahl 
(120101 ) for extending the extinction curve to the visible (orange dots). For comparison wc add a value for Ay derived from the X-ray 
spectrum of Sgr A* (Porquet et al. 2008) (viol et pent agon). In t he sam e way wc use the gas in front of Sgr A East jSakano ct al. 200^) 
(pink cross). The extinction curve differs from lRieke fc Lcbofs kvl II1985I): iRiekc ct al. (1989) (blue triangles), the results of which arc only 
partly based on GC data. The GC data differs also from the ICardelli et al.l ( 119891 ) curve (fitted to =3.1 obse rvations; green dashed 
line). The optical extinction towards the GC is uncertain. We use three curves of the type used bv ICardelli et aP lfl989i ) with Rv = 2.3, 
3.1, 3.9 (dotted gray lines from top down) for showing the possible range. 



IRawlings et all ()2003[ ) and lWhittet et al.l (119971 ) did not 
detect H2O ice in sight lines consisting of diffuse inter- 
stellar material with Ay w 10. Thus, these features 
are likely only visible in sight lines through molecu- 
lar clouds and not in purely diffuse ex tinction regions 
(jWhittet et all [l997l : iChiar et al.l[2000l ). This view is 
supported by t he CO clouds in the Galactic arms in front 
of the GC (S utton et al .l [l990l). This CO is also visible in 
the I SO-SWS spectrum ILutz et al.|[Tl96t iMoneti et al.l 
|2001|) . The detection of molecular cloud features towards 
the GC is not unique, ice features are also detected to- 
wards the Quintuplet cluster at a distance of about 12' 
to the GC, where they ar e slig htly weaker compared to 
the continuum extinction (jChiar et al.ll2000D . 

Other extinction features, lik e the strong aliphat ic hy- 
drocarbon feature at 3.4 ^m (jWillner et al.|[l979l ). are 
caused by diffuse du st and are thus a general visible 
featu re of extinction (jChiar et al.ll2000t IRawlings et al.l 
I2003f l. The feature at 3.4 /im is about a factor of two 



stron ger towards the GC than in the local diffuse extinc- 
tion ( Rawlings et"aII[2o5l ICao et al.l[20Tol ). The 3.4 ^m 
feature is typica l for diffuse dust. In molecular c louds it 
is not detected (jPendleton fc Allamandolall2002[ ). 

In order to better constrain the shape of the extinction 
curve, we use the continuum of the ISO-SWS spectrum, 
see Appendix [B] and Figure [H 

For the silicate feature at 9.7 /.tm we obtain an optical 
depth of Arsi9.7 — 3.84 ± 0.52 relative to the contin- 
uum at 7 /im from our interpolated extin ction curve. 
The d epth is similar to that o btained by IChiar et al.l 
(|2000l) of Arsis.T = 3.46 and by 'Roche & Aitken" fl985| 
of Tsi9.7 — 3.6. According to van Brecmcn ct al. (201(1) 
the shape of 9.7 /im silicate feature of the GC is iden- 
tical to sight lines with diffuse extinction and slightly 
different to sight lines with molecular clouds. Using 
our broadband extinction values (Appendix [C]) we ob- 
tain ATRi3.7/J^fJ - K) = 0.70 ± 0.10, consistent with 
(IRoche fc AitkenI [19851) . This is more than the value of 
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Arsi (). jlE(J — K) — 34 in nearly all other diffuse sigh t 
lines (|Roche fc AitkenI 119841: Ivan Breemen et alj [2010 ) . 
Sight lines with molecular clouds als o have rs, ■^.T/EjJ — 
K) = 0.34 or even smaller values (Ivan Breemen et all 
l2010t) . The large ATsig.7/E{J ~ K) towards the GC is 
probably caused by an abnormally high Arsi 9.7. A larger 
silicate dust to carbon dust ratio in t he inner Galaxy 
compared to the local diffuse m edium ([Roche fc AitkenI 
119851 : Ivan Breemen et all |201C1[) could explain the high 
ATsig.7/E{J—K) towards the GC. With a higher silicate 
abundance At^^^^^/ E{J — K) towards the GC should be 
identical to the local value. However, it is twice the local 
value. Porous dust grains ca use both a stron g 3.4 /itm 
and a strong 9.7 /zm feature ()Gao et al.l |2(310D. Poros- 
ity is also one element of the dust model of lZubko et al.l 
(|2004[ ). which best fits our continuum extinction data, 
see Section l5^ 

In the continuuni , CO2 is visible at 15 /zm 
(jGerakines et al.lll999f) . It is, however, much weaker than 
in the observ ations of e xtinct ion in molecular clouds ob- 
served by Kn ez et all ()2005[ ). In contrast, the silicate 
feature at 18 /im towards the GC is much stronger than 
t he CO2 feature, see Figur e El 

IRieke fc Leb"ofek^ (|l985l ) and IRosenthal et~all (pOOOl ) 
can fit all their extinction values up to 7.5 fiva. with a 
single power-law. This single power-law extinction is in- 
consistent with the flattening of the extinction curve in 
our data at around 4 /im, see Section 14.31 In the con- 
tinuum, there are no strong extinction features apparent 
between 3.7 and 8 /im, see Figure [51 Some weak fea- 
tures due to H2O, NH3, CO2, HCOOH and aliphatic hy - 
dro carbon are tentativ ely identified bv lLutz et al.l ([1996f ) 
and IChiar et al.l ([200C0 . All the features are weaker than 
r = 0.15. Thus, these features are too weak to explain 
a significant part of the extinction difference between 
the measured extinction and a single power-law. There 
is some indication for a steeper slope around 3.9 /zm 
an d a flatter one at longe r wavelengths. Also, the data 
of iNis hivama e t all ([2009D and most theoretical models 
('Section l5.8p exhibit a smooth flattening in the transition 
region. 
Zasowski et al.l 



TABLE 6 

Relative broadband extinction values derived from the 
GC extinction curve 



20061): IGao et al.l 



(20091: 
( 20091) : 

(EootI) 



IJiang et al.l 
iFlahertv et al.l 
INishivama et all 




Ronian -Zunig a et al.l .. .. , , , . . 

and [indcbctou wet al.l ([2005D also detected a flattening 
of the extinction at MIR wavelengths towards different 
reg ions of the Galax y . Som e of t hese observations, suc h 
as, INishivama et al.l ([20091 ) and IZasowski etHI ([20090 . 
targeted diffuse extinction in the bulge and the Galactic 
disk. Therefore, the flattening of the curve in the MIR 
is not caused by the molecular cloud in front of the GC. 

The extinction curves of other galaxies in the MIR 
are more difficult to interpret, because the extinction 
in other galaxies is not likely to be caused by an uni- 
form foreground screen. Howe ver, since a mixed model 
([Forster Schreiber et al.|[2"001|) produces a flattening of 
the extinction curve at higher extinction, while in the 
GC a flattening at smaller extinction is observed, it is 
still possible to find qualitative signs for the GC extinc- 
tion curve in the spectra of other galaxies. The fact that 
a single power-law extinction is slightl y preferred when 
fitting the continuum in some ULIRGS ([Tran et al.ll200lD 
could also partially be caused by a GC-like extinction 
and the mixed model. On the other hand, iThuan et al.l 
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VIRCAM Y 


4.634 ± 0.103 


NACO J 


3.051 ± 0.069 


NACO H 


1.737 ± 0.027 


NACO Ks 


1.000 ± 0.000 


NACO L' 


0.450 ± 0.053 


NACO M' 


0.391 ± 0.094 


IRAC 1 


0.547 ± 0.052 


IRAC 2 


0.396 ± 0.082 


IRAC 3 


0.340 ± 0.094 


IRAC 4 


0.383 ± 0.116 



Note. — The extinction ratios arc obtained from the average of two 
calculations for the effective filter extinction as described in Appendix 
[U] using Vega-like stars with A^r-y — 1 and Aur-y — 5. The errors 
are made up of half of the extinction ratio difference between the two 
different Braekett-7 extinction and the extinction curve errors. We as- 
sume that the NIR extinction follows a power-law of a — — 2.11±0.06. 
Broadband filters of other instruments can have another transmission 
curve. This can result in an additional error of up 3.5%. For further 
filters, see Appendix fC] 



(fT999l) found an extinction curve similar to our results for 
the extremely metal poor galaxy SBS 0335-052, also us- 
ing continuum emission. In the case of the central region 
of M82, Forster Schreib er et"all (|200H ) found, using hy- 
drogen lines, that an extinction curve similar to the one 
in the GC provides a better fit to the data than a single 
power-law. All in all, it is likely that the GC extinction 
([Lutz et al]|1996t ). considered to be unusual at the time 
of discovery, is in general the more widely spread type of 
extinction. 

5.7. Broadband extinction curve 

Since the GC extinction is consistent with many ex- 
tinction measurements outside of the GC (Section 15.4! 
and 15. 6|), the extinction curve towards the GC (Ap- 
pendix [B]) is also useful for other sight lines. We derive 
from this extinction curve broadband extinction values 
(Appendix [Cj . In principle it is necessary to know the 
extinction amplitude and the intrinsic object spectrum 
for deriving precise broadband extinctions. The error 
caused by neglecting these aspects is below 2 %. On the 
other hand, the error due to transmission differences be- 
tween slightly different, but equally named broadband 
filters can be up to 2.5 % (Appendix ICl) . We give in Ta- 
ble [H broadband extinctions ratios relative to the NACO 
Ks-band. For extinction values of many additional filters 
and for different source spectra, see Appendix [O 

5.8. Dust models 

Extinction curves, together with the emission spec- 
trum of dust, as well as the elemental abundances and 
depletions, constrain the properties of interstellar dust. 
Ideally, the dust model should also be plausible with 
regards to the formatio n and destruction of the dust 
(I Compiegne eFal] [20TTI ). We use published extinction 
models which fulfilled most constraints at the time of 
their publication and compare them with the extinction 
curve towards the GC. 

The classical grain model ([Mathis et al.lll977l) is com- 
posed of silicate and graphite grains, where both follow 
a power-law size distribution with a lower and an up- 
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Fig. 9. — Comparison of our data with other data and mod- 
els. To our extinction measurements derived from hydrog en lines 
(red b oxes) we add the inner bulge observations of Nishiva ma et al.l 
1120091 ) (blue circles) scaled to our ABr7. The models (lines) 
are from Wcingartncr & Drainc (2001); Zubko ct al. ( 2004): iDweH 
([2004) and Voshchinnikov ct al. (2006). Everything apart from the 
lowest data is shifted for better visualization. 



per cut off for the ex c lusion of very large and very small 
grains. iLi fc: Draind ()2001|) improved this simple model 
to to size distributions which include even smaller grains 
in order to account for the PAH emissions by interstellar 
clouds. We use a slight variant o f this dust model consid- 
ered by iWeingartner fc Draind ()2001[ ) (Wg- model) , who 
also consider different Ry- This model consists of a tri- 
modal carbonaceous grain size distribution and a simple 
silicate grain size distribution. With increasing Ry, the 
number of small grains decreases while the maximum 
grain size increases. 

Similarly to most later models the Wg-models do not 
have the H2O features at 3 /im and 6 /im. For constrain- 
ing the dust properties, independent of the features, we 
exclude the four lines concerned in our quantitative com- 
parisons. We also exclude the lines beyond 8 /im because 
the silicate feature is not modeled in all works and in 
no case matches our silicate extinction towards the GC. 
In addition, excluding these features we are also more 
likely to match other extinction measurements in the 
infrared, such as , for e xample [indcbctouw et al. (2005); 
INishivama et al.l (j2009D . where the features could not be 
measured due to too low spectral resolution. We fit all 
models to the data by adjusting the extinction curve with 
a global scaling factor. 

All Wg-models have a relatively flat NIR slope a > 
— 1.7 and hence fit the data badly in the NIR (A < 
2.8 /im), see Table H They differ in the MIR from each 
other. There the models with higher Ry are fiatter. 
On the one hand, when fitting all lines the models with 
Ry > 4.0 have a worse than the Ry =3.1 model 
because of the fiatter slope of the Ry > 4.0 models in 
the K-band, which is due to the smaller extinction errors 



there more important than the MIR. On the other hand, 
the Ry = 3.1 model i s incompatible with the data of 
INishivama et al ' (2009^ in the MIR, se e Figure H Com- 
pared with the Nishiva ma et al.l (j2009t ) data, the models 
with Ry > 4.0 are a better fit, although none of the Wg- 
models show such a strong and sharp fiattening from NIR 
to MIR as observed. 

One possibility for more comple x dust models is the 
inclu sion of voids in the dust grains (jVoshchinnikov et al.l 
120061) . In this model, mainly two size distributions of 
porous silicate grains and one of small graphite grains 
contribute to the extinction. Still, the extinction features 
shortward of 7 /tm are unexplained. This model has the 
advantage that the NIR slope between J and K depends 
on the porosity. None of their three models with different 
porosity and also no linear combination of these models 
can fit our data, sec Figure [5] and Table [71 

IZ ubko et al. ^ ( 2004, ) obtained dust models by fitting 
grain size distribution to an Ry = 3.1 extinction curve, 
dust emission and elemental abundances. They use the 
extinction dispersion between different sight lines as their 
error estimate. Thus, the extinction has a relatively 
small weight which could be part of the reason that var - 
ious dust models can fit their data. IZubko et all (|2004[ ) 
used different element abundances which have some im- 
pact on the result in the case of carbon abundances. We 
only use their solar abundances. More likely, the real 
solar abundances (Asplund ct al. 2009) are smaller than 
those used by (Zu bko et al.i 2004). However the abun- 
dances in HII regions are h igher than the solar abun- 
dances (|Asplund et al.ll2009() and at least the B-star car- 
bon abundance is too small fo r causing the observed mag- 
nitude of extinction (Li 20 0^1). 



The simplest modef of Zubko et al.l ([2004!), BARE- 
GR, consists of only graphite, PAHs and silicates grains 
and as such is similar to the Wg-models. BARE-GR 
only uses graphite grains smaller than 300 nm while the 
Wg-models use grains up to 0.8 //m and 5 fim respec- 
tively. Since grain sizes around 300 nm cause the steep- 
est infrared slopes in contrast to grains of other sizes 
(jMoore et al.l l2005: Gao et al.ll200 9V the NIR extinction 
slope of the BARE-GR model is steeper than that of 
Wg-models. OveraU, the BARE-GR model can fit the 
data with /d.o.f. = 18.9/14 relatively well. But this 
model does not have the H2O features in contrast to the 
ISO-SWS spectrum. Because this model nearly does not 
flatten at all in the MIR, it c annot fit the inner bulge 



data of INishivama et al 



The COMP models of 



(l2009l). see Figure la 



Zubko et all (|2004D use, in addi- 



tion to silicate and carbonaceous grains, composite grains 
consisting of silicates, organic refractory material, water 
ice and voids. These models seem to be more promis- 
ing for the GC, because they contain water ice which is 
visible in extinction features towards the GC. The best 
fitting subtype is AC in which the carbon is mainly amor- 
phous. The composite particles are, with sizes up to 
0.8 /im, bigger than the other particles. The other sub- 
types of composite models have a steeper slope, at least 
in the MIR, and hence make a worse fit to the data. The 
COMP-AC-S model has a x^/d-o.f. = 6.4/12, slightly 
better than the BARE-GR-S m odel. In contrast to the 
other models, it also matches the lNishivama et al.l (|2009f ) 
data, sec Figure [9l 
,Dwek (,2004 ) added metallic needles to the BARE-GR 
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TABLE 7 

Fit goodness of different dust models 



Author model total x^/d.o.f. best fit x'^ /points of NIR best fit /points of MIR 



WgOl Rv =3.1 114.3/14 106.6.3/6 7.7/9 

WgOl case B Rv =4.0 228/14 178/6 51/9 

Zu04 BARE-GR-S 18.9/14 4.4/6 14.5/9 

Zu04 COMP-AC-S 6.4/12 4.7/4 1.7/9 

Dw04 30.9/14 21.7/6 9.2/9 

Vo06 p= 0.18 111/14 66/6 45/9 

Vo06 p= 0.6 117/14 111/6 5.2/9 



Note. — We fit the different dust models to our data by scaling the full extinction curve. The NIR data are 
between 1 and 2.8 A^m, the MIR lines are between 3.7 and 7.5 fim. In x^/points we give the number of data 
points used for the f i tting which a r e with i n the NIR and MIR. The m odels are from IWeingartner Sc Draind 
ll2001l) : lDwei3 ll2004l) : [Zu"bko et al.l II2004I' ) : rVoshchinnikov et al.l i200B) . 



model of IZubko et al.l ()2004[ ). The extinction caused 
by the needles rises with wavelength up to 8 /xm, thus 
flattening the Wg-model between 3 and 8 /um, see Fig- 
ure [9l The extinction due to needles drops fast long- 
ward of 8 /im for avoiding a weakening of the silicate fea- 
ture. This model does not fit the data with x^/d-o.f. — 
30.9/14. The inconsistency is caused by an overly grad- 
ual transition to a flatter slope around 3 /im, while the 
observed slopes changes faster, see Figure 13 While it 
might be possible to fit our data better with some dust 
model variant which consists of metallic needles as well 
as simple carbonaceous and silicate grains, finding a fit 
to the data will not be easy, since the extinction due to 
both components varies so slowly with wavelength that it 
is difficult to obtain the observed sharp fiattening around 
3 fim. 

AU together, the COMP-AC-S model of IZubko et all 
(|2004l ) seems to be the best model for extinction towards 
the GC because of its best /d.o.f. and the presence of 
the H2O features. However, the lack of H2O ice features 
in m ost sight lines (jRawlings et al.l [20031 : TWhittet et al.l 
Il997f ). while MI R-flattenin g is ob served in m any regions 



Zasowski ot a l.ll2009 : 
2009; Fla hertv et al 



Jiang et al. 2003, 2006: lGao et alj 
2007: .Roman-Zufiiga et all 120071: 



Nishivam a et al.l l2009t llndebetouw et al.lT2005i r renders 
the idea that composite particles glued by ices are re- 
sponsible the fiat MIR extinction questionable, also for 
the GC. This is because if composite particles can exist 
only together with ices, the other regions should show a 
steeper MIR extinction than the GC. However, because 
the sight lines with extinction probed in the MIR do not 
overlap with the sight lines which are tested for H2O ice, 
it is still possible that the ice feature could exist in most 
sight lines with flatter MIR extinction. 

An alternative model could be that something else 
aside from ices produces the fiat extinction in the MIR 
towards the GC and elsewhere, and towards the GC ad- 
ditional pure ice grains produce the extinction features. 
However, all other models tested here have the problem 
that they do not have at once a steep NIR extinction and 
a fiat MIR extinction. It is perhaps possible to change 
the shape of the extinction curves within the models, 
without changing the type of particles in the model. 
However, this seems difficult for the pure carbonaceous 
and silicate grain model, because a smooth particle size 
distribution of carbonaceous and silicate grains should 
also produce a smooth extinction curve. It might be pos- 



sible, on the other hand, to change the extinction shape 
only around 3 /.tm by omitting only H2O ice from the 
composite particles in the COMP-AC-S model. If this 
does not change the other extinction properties of the 
dust such model could fit the extinction also in other 
sight lines. 

For solving the issue, further modeling seems to be 
necessary, maybe concentrating first on the detailed ex- 
tinction curve towards the Galactic Center. In the fu- 
ture, similarly detailed extinction curves in other sight 
lines would also be very useful. Especially promising 
could be searching for the H2O ice feature in sight lines 
which are already probed by broadband measurements. 
Also, adding more extinction measurements in between 
for testing the commonness of the sharp transition to a 
fiatter extinction in the MIR might be very interesting 
for further study. 

6. SUMMARY 

The simplest way to derive an extinction curve is to ob- 
serve a well understood object at different wavelengths, 
including at least one extinction-free wavelength. The 
minispiral in the GC fulfills these conditions. It is a HII 
region and Case B is valid. As extinction-free emission we 
use the 2 cm radio continuum observed with the VLA. In 
the infrared, we obtain line fiuxes between 1.28 /xm and 
2.17 /im from SINFONI. For lines at longer wavelengths 
of up to 19 /im we use ISO-SWS observations. We obtain 
the following results: 

• By interpolation we obtain A2.i66pm = 2.62 ± 0.11 
as the average 2.166 /im extinction of the ISO field 
of 14" x20" about Sgr A*. 

• Using the extinction map of lSchodel et"all ([20lol) 

for the relative spatial extinction we obtain, for the 
direction towards Sgr A*: Ar = 4.21±0.10, Aks = 
2.42 ± 0.10 and Al- = 1.09 ± 0.13. 



' iSchodel et al.l ()2010[ ) measured the total luminos- 
ity modulus of the red clump consisting of the ex- 
tinction and distance modulus towards the Galac- 
tic Center. Since we measure the extinction in- 
dependent of the distance, the combination of 
ISchodel et al.l ()2010[ ) with our extinction yields the 
distance Rq, the distance to the Galactic Center: 
we obtain Rq = 7.94 ± 0.65 kpc in agreement with 
current measurements. 
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• The extinction in the NIR (1.2 to 2.8 /im) is well fit- 
ted by a power-law of slope a = —2.11 ±0.06. This 
law is steeper than the value of about a ss —1.75, 
which was mostly report ed in the lit erature be- 
fore 2004, see for example (jDraineiri9890 . However, 
since 2005 most publications about diffuse extinc- 
tion t owards the bulg e and Galactic disk (for ex- 
ample iStead fc Hoard (|2009[ )) yield a steeper law. 
We obtain a = — 2.07±0.16 as weighted average of 
all publications about diffuse extinction since 2005. 

• At longer wavelengths, several extinction features 
are visible, such as H2O at 3.1 /im and silicates 



at 9.7 /im. Even aside from such features, the ex- 
tinction is higher than expected from extrapolating 
the NIR power-law. Our data agree well with sev- 
eral IRAQ publications for the inner bulge, such as 
iNishivama et al.l ([2009) for example. Because we 
are able to use many lines, it is apparent that the 
change in slope is sharper and stronger than in any 
of the extinction curves pro duced by pure carbona- 
ceous and silicate grains (jWeingartner fc Draind 
I2OOII ). The best fitting model adds composite 
parti cles which contain also H2O ice (jZubko et aP 
[2OOI . 
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APPENDIX 



A: CORRECTING FOR THE EXTINCTION LAW FLATTENING BIAS DUE TO INHOMOGENEOUS EXTINCTION 
Here, we estimate the extinction law flattening bias due to inhomogeneous extinction. 

We make the assumption that the true extinction map AjjX), f or rcsohition element i and wavelength A, has the 
same relative spatial distribution of extinction as the Scho deTet al.l ([2010) map, such that we can obtain the extinction 
map at a given wavelength A via a simple scaling of an (as yet unknown) factor x(A): 



A,{X) ^ x{X) X Aisch (Al) 

We t hen derive an unextincted Brackett-7 flux map from our observed Brackett-7 flux map, using the extinction 
map of iSchodel et al.l (|2010[ ) (scaled to our measured Brackett-7 extinction) as model for the spatial inhomogeneity in 
the extinction at this wavelength. For resolution element i, the unextincted flux map can be written as: 

F,{X^ Br7)u„ext = F,{X^ BtjU, * ^qOA.„ Aim.p),^,_ ^^2) 

Here we use as scaling factor, xq, our integrated n ieasurement o f the B rackett gamma extinction (Section l4.2[) . divided 
by the (observed) flux-weighted extinction of the ISchodel et al.l (|2010f ) map: 

A f Y^ieiso^ij^ = Br7)cxt A.sch ^ ^ , , 

2^0 — /i(A — or 7jiiitegrated, measured I tp / \ 5 \ I • V^'^J 

V Z^ielSO-^H^ = jext / 

We use the unextincted Brackett-7 flux map as a model for the spatial distribution of the intrinsic (unextincted) flux 
at all wavelengths (i.e. the correct distribution of relative flux weights). 
The true (extinction-law-conserving) integrated extinction should be flux-weighted by the unextincted flux, not the 
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Fig. 10. — Measured and unextinced MIR eontinuum towards the GC. The black line is the observed ISO-SWS spectrum. We unextinct 
the continuum around the hydrogen lines by the hydrogen Une extinctions there (red boxes). We interpolate these points to obtain the 
unextincted continuum (blue line). 



observed flux: 

Our extinction measurement, on the other hand, measures rather the ratio of integrated fluxes: 

^(A)intogratod,mcasurcd = "2.5 X log „ , ^ 5 ^ (A5) 

VZ^^GISO A = -t!r7ju„oxt/ 

In this equation, areas with smaller extinction have higher observed fluxes relative to the unextincted flux and thus 
are given a higher weight. Upon integration this then leads to a smaller measured extinction than the true integrated 
extinction of Equation IA4I 

To estimate this bias, we simulate our extinction measurements, weighting the intrinsic extinction, x{X)AiSch, by 
the unextincted flux weights derived from the Brackett-7 map: 

A(A)i„tegrated.simulated = "2.5 X log ^ilSlSO ^ (^g) 

V EielSO^ii-^ = Br7)uncxt / 



For every SINFONI line at wavelength A we find an a;(A) which, when substituted in Equation IA6[ results in our 
measured extinction A(A)intcgrated, measured- The calculation is carried out over the same area (ISO) used in the 
extinction measurement for each line. The bias at a given wavelength A is then: 

AA{X) = A(A)i„teg, ■ated.true -^('^) integrated, simulated 

(A7) 

This correction must be applied to A(A)intogratod, measured in order to obtain the correct (extinction-law-conserving) 
extinction. The biggest correction occurs for the highest extinction (Paschen-/?). Even for this line the correction is 
only 0.08 mag, smaller than the extinction error of 0.11 mag. 

B: INTERPOLATING THE EXTINCTION CURVE 

We use the continuum of the ISO-SWS spectrum in order to constrain the shape of the extinction curve in the 
MIR. The unextincted continuum cannot be modeled with a single simple model such as a single black- or graybody. 
Instead, we correct the measured continuum close to the hydrogen lines by using the extinction of the lines. We then 
interpolate the resulting points using simple models for sections of the spectrum, described in the paragraph below. 
After this, we obtain the extinction curve by dividing the interpolated extinction-corrected continuum by the observed 
continuum. By applying this method, we use the extinction of the hydrogen lines to obtain extinction values in regions 
where no hydrogen lines are available. 

Between 2.62 and 7.6 fim we fit the extinction corrected continuum points with a second order polynomial in log-log 
space. At longer wavelengths however the interpolation is less certain, due to t he fact that there are fewer lines, larger 
extinction errors and silicate emission in the GC region (jKemper et al.l [20041 ). For the continuum up to 12 fim we 
use the GC continuum of .Kemper et al.. (,2004) . which contains silicate emission, scaled such that it fits our two data 
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Interpolated infrared extinction curve 








UrS. \Kj\^ ) 


1.00009 


13.1075 


0.3 


1.00984 


12.8529 


0.2925 


1.0197 


12.6033 


0.285 


1.02964 


12.3586 


0.2776 


1.03969 


12.1186 


0.2701 


1.04983 


11.8832 


0.2626 


1.06007 


11.6525 


0.2551 


1.07041 


11.4262 


0.2477 


i.UoUoD 


1 1 ori/i Q 
ii.zU4o 




1.0914 


10.9867 


0.2327 


Note. 


— Our extinction curve. The second column is the average extinction 


towards the central 


14" x20" 


of the GC. The errors includes all effects. Some effects matter only for 


comparison over big 



wavelength ranges. Locally the error is smaller for most of ISO-SWS range (2.6 to 27 fim). This curve 
can be used for A(X) in Formula IC II for obtaining the extinction for other objects and filters. Scaled to 
other absolute extinctions it can be also useful outside the Galactic Center. The full table is available 
in the electronic edition. 



points in this wavelength range. A 240 K blackbody is a good approximation to the extincted continuum between 12 
and 15 fim. In this range there are no significant extinction features. We extrapolate this blackbody up to 26 /im, 
scaling it such that the unextincted continuum close to both lines in this range is fit by the blackbody, and use it as 
a model for the continuum between 12 and 26 /im. We combine the three different continua with smooth transitions 
around 7 //m and 12 /xm, see Figure [TUl 

In a similar way, we make a smooth transition to the a = —2.11 power-law shortward of 2.75 /im. We continue 
the interpolation by a a = —2.02 power-law shortward of 1.28 /xm We also correct unphysical jumps in the extinction 
curve, such as the one at 4.07 /im, due to calibration errors of the ISO-SWS data. For this we adjust the extinction 
curve locally with linear interpolation. We interpolate the extinction curve we obtain, excluding spectral lines, with 
splines (Table [5]). The spectral resolution of this interpolation is not high enough to fully resolve all features, such as, 
e.g., the CO feature at 4.7 /im. 

C: EFFECTIVE BROADBAND FILTER EXTINCTIONS 

A complication in correcting for extinction in broadband flux measurements is that the extinction that should be 
associated with a given filter varies with both the intrinsic spectrum of the source and the strength of the extinction. 

Since we use spectral lines in this paper to constrain the extinction law, the wavelengths to associate with the 
extinction measurements are well-defined. Also the interpolated extinction (Appendix[B]) has well defined wavel engths. 
Thus, we can explicitly calculate the extinction values for a given broadband filter, like lEspinoza et all ()2009[ ). using 
the equation: 

AHiter = -2.5 X log J ^ F,iX) S (X) dX ) ■ 

Here Fa(A) is the intrinsic fiux density of the source, which when multiplied with the wavelength. A, is proportional 
to the number distribution of photons with wavelength. S{X) is the wavelength dependent throughput of telescope, 
instrument and atmosphere. We also define here R{X) = lO^"'' '*^^'^)) as a "reddening factor" to separate the extinction 
from the intrinsic fiux. 

We use our interpolated extinction curve for A{X) (see Sect ion [Bland Table[5]) in Formula lCll to calculate the effective 
broadband filter extinction values for a range of commonly used infrared filters. The results are presented in Table O 
As source spectrum we use primarily a blackbody of 9480 K. For stars of this temperature, there is no major difference 
between the real spectrum and the correct blackbody in the infrared. The use of such a star also has the advantage that 
a similar star is used for the definition of the Vega magnitude systern . In addition, we also compute the effective filter 
extinctions for a bright Sgr A* like spectrum (lEisenhauer et al1l2005t Gillessen et al.|[2C)06D which follows a power-law 



with a slope /3 = 0.5 in i^Li, (jHornstein et al.| [2007: Dodds- Eden et alll2009D as example of a very red source. For S{X) 



we use the atmospheric transmission multiplied with the instrument transmission. For all Paranal filters we use the 
atran transmission models (Lord 1992) from Cerro Panchon for the atmosphere, with airmass 1 and 2.3 mm water 
vapor column from the Gemini web sit^l- For the VISIR filters we use the filter transmissions from the instrument web 
site0. We obtain the VIRCAM filters from the instrument web sit43- For the NIRC2 filters, which we obtain from 
the instrument web siteF^. we use a Mauna Kea atmosphere of airmass 1.5 and 1.6 mm water vapor column from the 
Gemini web site (see footnote [S]) . For these instruments only in case of VIRCAM the wavelength dependent quantum 
efficiency and mirror refiectivity are available. We use them for our calculation. However, including QE and mirror 

http.//vfwvf.gcrnini.edu/ . q--nodc/10789 l";" http://www.eso.org/sci/facilities/paranal/instruments/vircam/inst/ 

http://www.cso.org/sci/facilitics/paranal/mstruments/visir/mst/mde>i.ii1g^^p^//^^^2.keck.hawaii.edu/inst/nirc2/mters 
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TABLE 9 

Extinction of broadband filters for different Brackett-7 extinctions 
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Ks 
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2.42 ±0.10 


NACO 


L' 


1.20 ±0.14 


0.23 ±0.03 


2.73 ±0.32 


1.09 ±0.13 


NACO 


M' 


1.05 ±0.25 


0.20 ±0.05 


2.37 ±0.57 


0.95 ±0.23 


VISIR 


PAHl (8.6 /im) 


1.73 ±0.54 


0.34 ±0.10 


3.83 ± 1.23 


1.60 ±0.50 


VISIR 


PAH2_2 (11.88 Atm) 


1.61 ±0.36 


0.31 ±0.07 


3.66 ±0.83 


1.47 ±0.33 


NICMOS 


llOM 


10.42 ± 0.19 


2.07 ±0.04 


22.86 ±0.40 


9.49 ±0.18 


NICMOS 


145M 


6.00 ± 0.11 


1.16 ± 0.02 


13.44 ± 0.25 


5.47 ± 0.10 


NICMOS 


160W 


4.89 ±0.11 


0.97 ±0.02 


10.61 ±0.25 


4.40 ±0.10 


NICMOS 


170M 


4.33 ±0.11 


0.83 ±0.02 


9.80 ±0.25 


3.95 ±0.10 


NICMOS 


222M 


2.50 ±0.11 


0.48 ± 0.02 


5.70 ±0.26 


2.28 ±0.10 


IRAC 


Band 1 


1.47 ±0.14 


0.28 ±0.03 


3.27 ±0.32 


1.31 ±0.13 


IRAC 


Band 2 


1.06 ±0.22 


0.20 ±0.04 


2.40 ±0.50 


0.97 ±0.20 


IRAC 


Band 3 


0.91 ±0.25 


0.17 ±0.05 


2.08 ±0.57 


0.83 ±0.23 


IRAC 


Band 4 


1.02 ±0.31 


0.21 ±0.06 


2.19 ±0.67 


1.01 ±0.30 


NIRC2 


H 


4.75 ±0.11 


0.93 ±0.02 


10.56 ±0.25 


4.31 ±0.10 


NIRC2 


K' 


2.73 ±0.11 


0.53 ±0.02 


6.15 ±0.26 


2.48 ±0.10 


NIRC2 


L' 


1.24 ±0.14 


0.24 ±0.03 


2.80 ±0.32 


1.12 ±0.13 


NIRC2 


Ms 


1.17 ±0.25 


0.23 ±0.05 


2.66 ±0.57 


1.07 ±0.23 


VIRCAM 


Y 


12.40 ±0.27 


2.40 ± 0.05 


27.77 ±0.58 


11.33 ±0.25 


VIRCAM 


J 


8.21 ±0.12 


1.60 ±0.03 


18.30 ±0.26 


7.48 ±0.11 


VIRCAM 


H 


4.68 ±0.11 


0.91 ±0.02 


10.41 ±0.25 


4.24 ±0.10 


VIRCAM 


Ks 


2.67 ±0.11 


0.51 ±0.02 


6.04 ±0.26 


2.42 ±0.10 


2MASS 


J 


8.26 ±0.12 


1.64 ±0.03 


18.14 ±0.26 


7.51 ±0.11 


2MASS 


H 


4.65 ±0.11 


0.90 ±0.02 


10.37 ±0.25 


4.22 ±0.10 


2MASS 


Ks 


2.58 ±0.11 


0.50 ±0.02 


5.84 ±0.26 


2.34 ±0.10 



Note. — Wc calculate filter extinctions for different Brackett-7 (2.166 /^m) extinctions. We use blackbodics and 
power-law source spectra. The error is the uncertainty due to the extinction error. The first column is the average 
extinction of stars in the central 14" x20" of the GC. The last column is the extinction towards the power-law source 
Sgr A* (i.e. using stars in the central r< 0.5" of the GC for scaling of Asr-f )■ 



reflectivity change the extinction only by less than 0.07% relatively. For the other instruments we assume that the 
throughput of the instrument is apart of the filter not wavelength d epe ndent within a filter. For 2MASS we use the 
full transmissions including the atmosphere from the project web sitem. We obtain the NICMOS transmissions from 
the instrument web siteHEI as in case of the IRAC transmission^. 

Although, strictly speaking, it is necessary to calculate the effective broadband extinction on a per spectrum basis, 
the differences in the effective extinction for different source SEDs are mostly relatively small, even for the high 
extinction of the GC. For example, the difference in effective extinction between a blue spectrum (Vega) and a red 
spectrum (Sgr A) is only 0.056 mag in the H-band and 0.026 mag in the Ks-band. 

We test the influence of difference magnitudes of extinctions on broad band extinctions. The H-band extinction 
derived from Asr-y = 2.62 is 4.65 mag, while the H-band extinction derived from Asr-i = 0.5 is 0.91 mag. Assuming 
linear scaling of the H-band extinction with the Brackett-7 extinction the expected H-band extinction (scaled up from 
Asr^y = 0.5) is 0.91mag x 2.62/0.5 = 4.77 mag. Thus, the deviation from linearity (the non-linearity) is 0.12 mag in 
this case. In the Ks-band the non-linearity is, with 0.029 mag, much smaller for the same value of Abti- We find that 
the non-linearity is of the order of the measurement error for Aband ^ 4. 

For most objects, for which the extinction is not significantly higher than in the GC, it is sufficient to use the 
extinction closest to the measured extinction in Table[9]renormalized to the correct absolute extinction. For very high 
extinctions or different filters it is necessary to calculate the filter extinctions from Formula ICll using the interpolated 
extinction curve given in Table |H1 

D: COMPARISON WITH OTHER METHODS FOR THE DETERMINATION OF THE NIR EXTINCTION SLOPE 

We derive the NIR extinction slope a from absolute extinction values at known wavelength. Another way commonly 
used to compute the near infrared extinction slope a is to use stellar colors, using an assumed wavelength to associate 
with the broadband filters (we hereafter call this the effective wavelength method). For exa mple, the following equation 
is commonly used to compute the extinction slope a from JHK(s) colors (jStead fc Hoarell2009.) : 



Ej-H I A 



^) -1 



Eh-Ks I _ 



Aff 



(Dl) 



http://www.ipac. caltech. edu/2mass/releases/second/doc/soc3_lbl. html#sl6 / / ■, ij. i, j /■ / vu ^- c-\ i ^ i 

,0 /, , /, ; /, In, nttp://ssc. spitzer.caltecn.edu/irac/calibrationmcs/spectralresponse 

http:/ /www.stsci.edu/hst/mcmos/design/nlters 
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where Aj is some assumed wavelength of the J filter, and so on for the other filters. Obviously, it is necessary to know 
the wavelength that can correctly be associated with the filters for the given source and exti nction in order to derive 
a. One possibility would be to use the isophotal wavelength (Aiso) (|Tokunaga fc Vaccall2005D : 

/F,(A)5(A)rfA 
^'^^-^ - JS{X)dX 

This means that Ajso is the wavelength at which the monochromatic flux F\( X;i^„) equals the mean flux in the passband. 
More commonly used for extinction purposes is the effective wavelength (jTokunaga fc Vaccall2005j ): 

^ JX^F^iX)S{X) dX 
Aoff uncxt JXFx{X)SiX) dX ^ ' 

The effective wavelength defined this way is the average wavelength of received photons, weighted by the number 
distribution of received photons at the detector, appropriate for photon counting detectors. This formula is used in 

two variants: for the calculation of X'fi the source spectra are not extincted. 

For the calculation of the other variant (A^gg^^) the source spectra are extincted (|Stead fc Hoarell2009l : iSchodel et al.l 
l2OT0h : 

^ /A^fA(A)^(A)i?(A)dA 
Acffcxt /A^^A(A)5(A)i?(A)dA ^ ' 

We test the accuracy of the methods concentrating on the two definitions of the effective wavelength. To do this 
we compare the explicitly calculated filter extinctions to those obtained with the effective wavelength method for an 
extincted, 9480K Vega-like blackbody. As test filters we use the NACO JHKs filters plus atmosphere and our a = —2.11 
extinction law. Using Formula ICll we obtain slightly higher extinctions for the explicit calculation than obtained by 
applying the extincted effective wavelength method (Formula lDip . with differences of ^/^s true — A' cffoxt = 0.05 for 
ABr7 = 2.5 (for lower extinction the difference is much less Axs tmc — ^ifsA'offoxt = 0.013 for ABr7 = 1). Similarly 
we obtain Aks true — ^ifs A' eff unext = —0.012 for ABr7 — 2.5 when we use the unextincted effective wavelengths. 

The reason for these discrepancies is illustrated via an exaggerated example in Figure [TTJ To obtain the correct 
value of the true extinction it is necessary to consider the photon distribution, integrated over the filter, of both the 
unextincted and of the extincted source. Yet the two effective wavelength methods each consider only one of the above 
photon distributions. In the case of the unextincted source: because a larger proportion of the photons is transmitted 
through the extinction at longer wavelength than at A[,g^j^j,^^, the true filter extinction is smaller than at X'^^^^^^^. 
The opposite is true at X'^g^^^. 

In summary both 'effective' wavelengths deviate from the true extinction, because these effective wavelengths are 
not the effective extinction wavelengths. As a result it is not correct to use these 'effective' wavelengths to calcu- 
late the extinction appropriate for broadband filters from our extinction curve. The true wavelength for extinction 
measurements Atrue is the wavelength a t wh ich the extincted object is extincted by the same amount of extinction 
as in the explicit calculation of Formula ICll Atrue de pends on the optical s y stem, object SEP, st rength of extinction 
and the shape of the extinction law. However, since iStead fc Hoa^ (|2009f ): ISchodel et al.l ()2010[ ). e.g., measured the 
broadband extinctions directly their extinctions are correct at least for the specific source, filter and extinction strength 
combinations used in the determination. 

More important is the issue of effective wavelengths for the determination of a fr om m easured broadband extinctions. 
To test this we calculate alpha, using either Aiso, A^ffm^^^^t -^offcxt For mula | D1 | fmeth o d of relative extinction). 
We also calculate a from the absolute Ah and Aks using A^^^^^, as is done in ISchodel et al.l (|2010f) . We then compare 
the obtained values of a with the input a = —2.11, see Figure [l31 

Since none of the methods obtains the true extinction values, it is not surprising that none of the methods obtains 
the true extinction slope a. The best method for obtaining a is using X'^ff^^^ and Formula ID II for extinctions ABr7 < 4 
there is less than 0.046 deviation of a from the input a = —2.11. The other methods can result in larger deviations of 
up to to 0.11 for the same range of extinction. 

The difference between the a obtained using X[^^ and A[,ff ^^^j^^ is only Aa «0.01. In contrast IStead fc Hoard ()2009l ) 
obtain a difference of around Aa =0.2. The reason for this could be the slightly different photometric system, but 
perhaps more likely, the difference between a 9480 K blackbody and a K2III spectrum. Given this big difference in 
Aa, it is unclear if using A^g u^ext best approximation to the true a for instruments and source spectra 

other than tested here. 
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Fig. 11. — Difference between the extinction at (extincted and unextincted) effective wavelengths, and the filter-integrated extinction. 
A source (blue Hne, showing unextincted flux) is observed with a broad filter (gray dots). Via an (unrealistic) extinction, the source is 
extincted to the orange line, fn this example the effective wavelength of the unextincted flux is in the part with high extinction, and the 
effective wavelength of the extincted flux is in the part with small extinction. The true reddening factor, however, is determined by the 
integrals over the filter range and has a value between the reddening factors at the two effective wavelengths. For more realistic extinction 
curves this effect is much smaller, but not always negligible. 




"Br/ 

Fig. 12. — Comparison of the true near infrared extinction slope a with the a obtained using common methods. We use the NACO 
JHKs fil ters with atmosphere and a 9480 K blackbody, simulating common measurements. For obtaining a f rom Ej_u/Eu_xs use 
Formula ID II The results are displayed as follows: (o rang e) the use of the isophotal wavelength in Formula IDll (green) the use of the 
effective wavelength of the unextincted flux in Formula ID II (black) the use of the effective wavelength of the extincted flux in Formula ID II 
and (blue) the computation of a from absolute Ajj and Axs , using the extincted effective wavelengths. The true slope is shown as the red 
line. 



